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ABSTRACT 

We  present  the  results  of  Spitzer  Infrared  Spectrograph  low-resolution  infrared  5-35  pm  spectroscopy  of  nearby 
ultraluminous  infrared  galaxies  (ULIRGs)  at  z  <  0.15.  We  focus  on  the  search  for  the  signatures  of  buried  active 
galactic  nuclei  (AGNs)  in  the  complete  sample  of  ULIRGs  classified  optically  as  non-Seyferts  (LINERs  or  H  n 
regions).  In  addition  to  polycyclic  aromatic  hydrocarbon  (PAH)  emission  features  at  6.2,  7.7,  and  11.3  p m,  the 
conventional  tool  of  starburst-AGN  separation,  we  use  the  optical  depths  of  the  9.7  and  1 8  pm  silicate  dust  absorption 
features  to  infer  the  geometry  of  energy  sources  and  dust  at  the  nuclei  of  these  ULIRGs,  namely,  whether  the  energy 
sources  are  spatially  well  mixed  with  dust  (a  normal  starburst)  or  are  more  centrally  concentrated  than  the  dust  (a 
buried  AGN).  Infrared  spectra  of  at  least  30%,  and  possibly  50%,  of  the  observed  optical  non-Seyfert  ULIRGs  are 
naturally  explained  by  emission  consisting  of  (1)  energetically  insignificant,  modestly  obscured  (Ay  <  20-30  mag) 
PAH-emitting  normal  starbursts  and  (2)  energetically  dominant,  highly  dust-obscured,  centrally  concentrated  energy 
sources  with  no  PAH  emission.  We  interpret  the  latter  component  as  a  buried  AGN.  The  fraction  of  ULIRGs  showing 
some  buried  AGN  signatures  is  higher  in  LINER  ULIRGs  than  in  H  n  region  ULIRGs.  Most  of  the  luminous  buried 
AGN  candidates  are  found  in  ULIRGs  with  cool  far-infrared  colors.  Where  the  absorption-corrected  intrinsic  AGN 
luminosities  are  derivable  with  little  uncertainty,  they  are  found  to  be  of  the  order  of  1012  L0,  accounting  for  the  bulk 
of  the  ULIRGs’  luminosities.  The  5-35  pm  spectroscopic  starburst/AGN  classifications  are  generally  consistent  with 
our  previous  classifications  based  on  3-4  pm  spectra  for  the  same  sample. 

Subject  headings:  galaxies:  active  —  galaxies:  ISM  —  galaxies:  nuclei  —  galaxies:  Seyfert  — 
galaxies:  starburst  —  infrared:  galaxies 


1.  INTRODUCTION 

Galaxies  that  radiate  very  large  luminosities  as  infrared  dust 
emission  (Zir  <  1012  Z0)  are  called  ultraluminous  infrared  gal¬ 
axies  (  ULIRGs;  Sanders  &  Mirabel  1996).  The  huge  infrared 
luminosities  mean  that  (1)  very  powerful  energy  sources,  star- 
bursts  and/or  active  galactic  nuclei  (AGNs),  are  present  hidden 
behind  dust;  (2)  energetic  radiation  from  the  energy  sources  is 
absorbed  by  the  surrounding  dust;  and  (3)  the  heated  dust  grains 
re-emit  this  energy  as  infrared  thermal  radiation.  In  the  local 
universe  (z  <  0.3),  the  ULIRG  population  dominates  the  bright 
end  of  the  luminosity  function  (Soifer  et  al.  1987).  The  impor¬ 
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tance  of  the  ULIRG  population  i.e.,  the  comoving  infrared  en¬ 
ergy  density,  increases  rapidly  with  increasing  redshift  (Le  Floc’h 
et  al.  2005).  Identifying  the  dust-obscured  energy  sources  of  nearby 
ULIRGs  is  thus  very  important,  not  only  to  unveil  the  nature  of 
nearby  ULIRGs,  but  also  to  clarify  the  obscured  AGN-starburst 
connections  in  the  high-redshift  (z  >  0.5)  ULIRG  population, 
where  detailed  investigations  are  difficult  with  existing  facilities. 

In  nearby  ULIRGs,  it  is  now  widely  accepted  that  the  domi¬ 
nant  energy  sources  are  compact  (<500  pc),  highly  dust-obscured 
nuclear  cores,  rather  than  extended  (>kpc),  weakly  obscured  star- 
burst  activity  in  the  host  galaxies  (Soifer  et  al.  2000;  Fischer 
2000).  The  most  important  issue  is  to  determine  whether  the  ob¬ 
scured  compact  cores  are  powered  by  very  compact  starbursts 
and/or  AGNs,  from  which  strong  emission  is  produced  in  very 
compact  regions  around  the  central  accreting  supermassive  black 
holes. 
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If  a  luminous  AGN  is  present  hidden  behind  dust  in  a  torus 
geometry,  as  inferred  for  the  classical  Seyfert  population  (e.g., 
Antonucci  1993),  a  large  amount  of  energetic  AGN  radiation 
can  escape  along  the  direction  perpendicular  to  the  torus  and 
form  the  so-called  narrow-line  regions  (NLRs;  Antonucci  1 993). 
NLRs  are  the  sources  of  strong  forbidden  emission  lines  seen 
in  the  optical-infrared,  and  their  relative  strengths  are  different 
from  those  of  a  starburst  because  the  spectrum  at  the  ionizing 
high-energy  UV  range  (A  <  912  A)  is  harder  in  an  AGN  than  in  a 
starburst.  An  AGN  obscured  by  a  torus-shaped  dust  distribution, 
with  a  well-developed  NLR,  is  thus  distinguishable  from  a  star- 
burst  galaxy  relatively  easily,  through  optical  spectroscopy  (classi¬ 
fied  as  Seyfert  2;  Veilleux  &  Osterbrock  1987)  or  high-resolution 
infrared  spectroscopy  (Genzel  et  al.  1998;  Armus  et  al.  2007). 

Since  the  gas  and  dust  in  an  AGN  have  angular  momentum 
with  respect  to  the  central  supermassive  black  hole,  an  axisym- 
metric  spatial  distribution  is  more  natural  than  spherical  geometry. 
In  this  case,  the  column  densities  can  be  high  in  certain  directions 
(the  torus  directions)  but  low  in  others  (along  the  torus  axis).  The 
classical  Seyfert-type  AGN  is  naturally  understood  as  the  popu¬ 
lation  where  most  of  the  torus  axis  direction  is  transparent  to  the 
AGN’s  ionizing  UV-soft-X-ray  radiation.  However,  compared  to 
the  classical  Seyfert  AGNs,  ULIRGs  are  known  to  contain  a  sub¬ 
stantially  larger  amount  of  gas  and  dust  in  their  nuclei  (Sanders  & 
Mirabel  1996).  For  a  fixed  angular  momentum,  the  increased 
amount  of  gas  and  dust  results  in  higher  column  densities  in  all 
directions.  If  high-density  gas/dust  blocks  the  bulk  of  the  AGN’s 
ionizing  radiation  at  very  small  radii  (<SC10  pc)  even  in  the  torus 
axis  direction,  then  no  significant  NLR  will  be  present.  It  is  very 
likely  that  any  luminous  AGNs  in  many  ULIRGs  exist  in  this 
buried  condition,  with  no  obvious  Seyfert  signatures  in  the  opti¬ 
cal  spectra  (  hereafter  “optical  non-Seyfert”).  Studying  only  AGNs 
with  well-developed  NLRs  provides  a  very  incomplete  picture  of 
the  AGN  population  and  will  miss  many  buried  AGNs,  whose  un¬ 
derstanding  may  be  essential  to  clarify  the  true  nature  of  ULIRGs. 

An  often-used  indicator  to  separate  a  buried  AGN  from  a  nor¬ 
mal  starburst  is  the  emission  from  polycyclic  aromatic  hydro¬ 
carbons  (PAHs),  seen  in  infrared  spectra  at  Arest  =  3-25  pm  in 
the  rest  frame  (Genzel  et  al.  1998;  Imanishi  &  Dudley  2000).  In 
a  normal  starburst,  UV-emitting  H  n  regions  and  molecular  gas 
and  dust  are  spatially  well  mixed.  PAHs  are  excited  by  far-U  V  pho¬ 
tons  from  stars,  and  strong  PAH  emission  is  produced  in  photo¬ 
dissociation  regions  (PDRs),  the  interfaces  between  H  n  regions 
and  molecular  gas  (Sellgren  1981),  unless  the  metallicity  is  very 
low  (Houck  et  al.  2004b).  In  a  normal  starburst,  the  equivalent 
width  of  the  PAH  emission  is  insensitive  to  level  of  dust  extinc¬ 
tion,  as  both  the  PAH  emission  and  the  nearby  continuum  are 
similarly  attenuated.  Thus,  a  normal  starburst  should  always  show 
large  equivalent  width  PAH  emission,  regardless  of  dust  extinc¬ 
tion.  Near  an  AGN,  the  PAHs  are  destroyed  by  strong  X-ray  ra¬ 
diation  from  the  AGN  (Voit  1992;  Siebenmorgen  et  al.  2004).  If 
the  PAHs  are  sufficiently  shielded  by  obscuration  from  the  X-ray 
emission  of  the  central  AGN,  the  PAHs  can  survive.  However,  the 
PAH-exciting  far-UV  emission  from  the  AGN  will  also  be  at¬ 
tenuated  by  the  obscuration,  so  in  a  pure  AGN  without  starburst 
activity,  the  PAH  emission  will  be  very  weak.  Instead,  a  PAH- 
free  continuum  from  hot,  submicron-sized  dust  grains  heated  by 
the  AGN  is  observed.  In  a  starburst/AGN  composite  galaxy,  PAH 
emission  from  the  starburst  regions  is  seen,  but  its  equivalent 
width  will  be  smaller  than  in  a  pure  starburst,  because  of  the  dilu¬ 
tion  by  PAH-free  continuum  produced  by  the  AGN.  Thus,  this 
PAH-based  method  can  be  used  to  disentangle  a  buried  AGN 
from  a  normal  starburst  at  an  obscured  ULIRG’s  core.  Since  the 
PAH  emission  features  are  spectrally  very  broad,  low-resolution  in¬ 


frared  spectroscopy  is  adequate.  In  addition,  ULIRGs  are  by  def¬ 
inition  bright  in  the  infrared,  so  that  this  infrared  PAH  technique 
can  be  applied  to  a  large  number  of  ULIRGs  in  a  reasonable  amount 
of  telescope  time,  and  thus  can  be  a  very  effective  tool  to  sys¬ 
tematically  identify  the  hidden  energy  sources  of  ULIRGs’  cores. 

Using  low-resolution  infrared  spectroscopy,  we  also  have  an¬ 
other  powerful  tool  to  distinguish  buried  AGNs  from  normal 
starbursts  in  the  cores  of  ULIRGs,  by  investigating  the  geometry 
between  energy  sources  and  the  gas  and  dust.  In  a  normal  star- 
burst,  the  energy  sources  (stars)  and  gas/dust  are  spatially  well 
mixed  (Fig.  la)  (Puxley  1991;  McLeod  et  al.  1993;  Forster 
Schreiber  et  al.  2001),  while  in  a  buried  AGN,  the  energy  source 
(the  central  accreting  supermassive  black  hole)  is  very  compact 
and  is  more  centrally  concentrated  than  the  surrounding  gas  and 
dust  (Fig.  lb)  (Soifer  et  al.  2000;  Imanishi  &  Maloney  2003; 
Siebenmorgen  et  al.  2004).  The  difference  in  geometry  is  re¬ 
flected  in  two  different  ways  in  the  observed  low-resolution  infra¬ 
red  spectra.  The  first  is  that,  while  the  absolute  optical  depths  of 
dust  absorption  features  in  the  3-10  pm  range  cannot  exceed  a  cer¬ 
tain  threshold  in  a  normal  starburst  with  mixed  dust/source  geom¬ 
etry  (Fig.  la),  they  can  be  arbitrarily  large  in  a  buried  AGN  (Fig.  lb) 
(Imanishi  &  Maloney  2003;  Imanishi  et  al.  2006a;  Levenson  et  al. 
2007).  The  second  is  that  a  strong  dust  temperature  gradient  is 
found  in  a  buried  AGN  (Fig.  2),  but  not  in  a  normal  starburst.  This 
dust  temperature  gradient  can  be  investigated  by  comparing  opti¬ 
cal  depths  of  dust  absorption  features  at  different  infrared  wave¬ 
lengths  ( Dudley  &  Wynn- Williams  1997;  Imanishi  2000;  Imanishi 
et  al.  2006a). 

In  this  paper  we  present  low-resolution  5-35  pm  spectra  of 
nearby  ULIRGs  at  z  <  0.15  to  search  for  the  signatures  of  op¬ 
tically  elusive  ( Maiolino  et  al.  2003)  buried  AGNs  in  the  cores  of 
ULIRGs  and  to  determine  the  energetic  contribution  of  buried 
AGNs.  Throughout  this  paper,  Ho  =  75  km  s_1  Mpc-1,  VLm  = 
0.3,  and  Da  =  0.7  are  adopted.  The  physical  scale  of  1"  is 
0.34  kpc  in  the  nearest  source  at  z  =  0.018, 1.72kpcatz  =  0.1, 
and  2.44  kpc  in  the  furthest  source  atz  =  0.15. 

2.  TARGETS 

Our  main  targets  are  optically  non-Seyfert  ULIRGs  at 
z  <  0.15  inthe/JMS  1  Jy  sample  (Kim  &  Sanders  1998).  In  this 
1  Jy  sample,  there  are  69  ULIRGs  at  z  <  0.15.  Based  on  the 
optical  classifications  by  Veilleux  et  al.  (1999a,  their  Table  2), 
28  ULIRGs  are  classified  optically  as  LINERs  and  20  ULIRGs  are 
classified  optically  as  H  n  regions.  These  48  (=28  +  20)  ULIRGs 
are  selected,  comprising  ~70%  of  the  1  Jy  sample  at  z  <  0.15. 
The  remaining  ~30%  ULIRGs  are  classified  optically  as  Seyferts 
(17  sources)  and  optically  unclassified  (4  sources;  Veilleux  et  al. 
1999a).  Table  1  summarizes  the  basic  information  and  IRAS- 
based  infrared  properties  of  the  selected  LINER  and  H  n  region 
ULIRGs.  Low-resolution  infrared  5-35  pm  spectra  for  this  com¬ 
plete  sample  are  presented  in  this  paper.  For  most  sources,  ancil¬ 
lary  ground-based  infrared  3-4  pm  spectra  are  available  (Imanishi 
et  al.  2006a). 

3.  OBSERVATIONS  AND  DATA  ANALYSIS 

Observations  of  all  the  48  ULIRGs  were  performed  using 
the  Infrared  Spectrograph  (IRS;  Houck  et  al.  2004a)  on  board 
Spitzer  (Werner  et  al.  2004).  All  four  modules,  Short-Low  2 
(SL2;  5. 2-7. 7  pm)  and  1  (SL1;  7.4-14.5  pm),  Long-Low  2 
(LL2;  14.0-21.3  /rm)and  1  (LL1;  1 9.5—38.0  pm)  were  used,  to 
obtain  full  5-35  pm  low-resolution  (R  ~  100)  spectra.  13  out 
of  the  48  sources  were  in  the  GTO  target  list  (PI:  J.  Houck, 
pid  =  105),  and  their  archival  spectra  were  analyzed.  The  re¬ 
maining  35  sources  were  observed  through  our  GO-1  program 
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PAH  +  continuum 
(weakly  obscured) 


PAH-free  continuum 
(highly  obscured) 


xtreme  starburst 
(Hll-regions  only) 


Fig.  1 . — (a)  Geometry  of  energy  sources  and  dust  in  a  normal  starburst.  The  open  star  symbols  indicate  “stars”  in  a  starburst.  The  energy  sources  (stars)  and  dust  are 
spatially  well  mixed.  Large  equivalent  width  PAH  emission  is  observed.  ( b )  Geometry  of  the  energy  source  and  dust  in  a  buried  AGN.  The  energy  source  (a  mass- 
accreting  supermassive  black  hole)  is  more  centrally  concentrated  than  the  surrounding  dust.  No  PAH  emission  is  found,  (c)  A  buried  AGN/starburst  composite.  The 
starburst  surrounds  the  central  buried  AGN.  The  observed  spectrum  is  a  superposition  of  PAH  +  continuum  emission  from  the  starburst  and  PAH-free  continuum  from 
the  AGN.  Since  the  buried  AGN  is  more  highly  obscured  than  the  surrounding  starbursts,  the  starburst  emission  generally  makes  a  strong  contribution  to  an  observed 
flux  even  if  the  starburst  is  energetically  insignificant,  (d)  A  normal  starburst  nucleus  (mixed  dust/source  geometry)  obscured  by  a  large  amount  of  foreground  dust  in  a 
host  galaxy.  This  geometry  can  produce  a  strong  dust  absorption  feature  whose  optical  depth  is  larger  than  the  maximum  threshold  obtained  in  the  mixed  dust/source 
geometry,  (e)  An  exceptionally  centrally  concentrated  extreme  starburst,  whose  emitting  volume  is  predominantly  occupied  by  H  n  regions.  Such  an  extreme  starburst 
produces  no  PAH  emission,  like  a  buried  AGN. 


(PI:  M.  Imanishi,  pid  =  2306)  and  another  GO-1  program  (PI: 
S.  Veilleux,  pid  =  3187).  The  GO-1  program  by  Dr.  Veilleux 
covered  only  SL2  and  SL1.  For  these  ULIRGs,  LL2  and  LL1 
spectra  were  taken  in  our  GO-1  program. 

For  our  program,  each  observation  consists  of  60  s  x  2  cycles 
for  SL  and  30  sx4  cycles  for  LL.  The  observation  was  per¬ 
formed  both  at  A  and  B  beam  positions,  and  so  total  exposure  time 
was  240  s  for  both  SL  and  LL.  For  the  programs  by  Houck  and 
Veilleux,  the  total  exposure  times  are  similar  or  smaller  for  bright 
sources.  Table  2  tabulates  the  observing  log  of  the  48  ULIRGs. 
The  slit  width  is  3.6"  or  2  pixels  for  SL2  (1.8"  pixel-1)  and  is  3.7" 
or  ^2  pixels  for  SL1  (1.8"  pixel-1).  For  LL2  and  LL1,  the  slit 
widths  are  10.5"  and  10.7",  respectively,  which  correspond  to 
~2  pixels  for  both  LL2  (5.1"  pixel-1)  and  LL1  (5.1"  pixel-1). 

The  latest  pipeline-processed  data  products  at  the  time  of  our 
data  analysis  (S 1 1  - 1 4,  pbcd  files)  were  used.  Frames  taken  at  the 
A  position  were  subtracted  from  those  taken  at  the  B  position,  to 
remove  background  emission,  mostly  the  zodiacal  light.  Then,  spec¬ 
tra  were  extracted  in  a  standard  manner  using  the  task  apall  in 
IRAF.3  Apertures  with  4-6  and  4-5  pixels  were  employed  for  SL 
and  LL  data,  respectively,  depending  on  the  spatial  extent  of  individ¬ 
ual  sources.  Then,  spectra  extracted  for  the  A  and  B  positions  were 
summed.  Wavelength  calibration  was  made  based  on  the  files  of 

3  IRAF  is  distributed  by  the  National  Optical  Astronomy  Observatories, 
which  are  operated  by  the  Association  of  Universities  for  Research  in  Astronomy, 
Inc.  (AURA),  under  cooperative  agreement  with  the  National  Science  Foundation. 


the  Spitzer  pipeline  processed  data,  named  “b0_wavsamp.tbl” 
and  “b2_wavsamp.tbF’  for  SL  and  LL,  respectively.  It  is  believed 
to  be  accurate  within  0. 1  pm.  A  small  level  of  possible  error  in 
wavelength  calibration  will  not  affect  our  main  discussions.  Since 
emission  from  all  ULIRGs  is  dominated  by  spatially  compact 
sources  at  the  observed  wavelength,  flux  calibration  was  per¬ 
formed  using  the  Spitzer  pipeline  processed  files  “b0_fluxcon 
.tbl”  (SL)  and  “b2_fluxcon.tbl”  (LL).  For  SL1  spectra,  data  at 
20bs  >  14.5  pm  in  the  observed  frame  are  invalid  (Infrared  Spec¬ 
trograph  Data  Handbook  ver.  1.0)  and  so  were  discarded.  For 
LL1  spectra,  we  used  only  data  at  20bs  <35  pm ,  because  the  scat¬ 
ter  is  large  at  A0bs  >35  pm.  No  defringing  was  attempted. 

For  flux  calibration,  we  adopted  the  values  of  the  pipeline 
processed  data.  We  made  no  attempt  to  recalibrate  our  spectra 
using  IRAS  measurements  at  12  and  25  pm,  primarily  because 
( 1 )  only  upper  limits  are  provided  for  IRAS  1 2  pm  photometry  in 
a  sizable  fraction  (>60%;  3 1/48)  of  the  observed  ULIRGs  and  for 
25  pm  data  in  some  fraction  (^20%;  10/48)  of  them  (Table  1), 
and  (2)  the  validity  of  recalibrating  IRS  fixed-aperture  slit  spec¬ 
tra,  using  the  large  aperture  photometry  of  IRAS,  is  not  obvious 
for  ULIRGs  that  possibly  show  a  non-negligible  fraction  of  spa¬ 
tially  extended  dim  emission  outside  the  IRS  slit.  Hence,  the 
absolute  flux  calibration  is  dependent  on  the  accuracy  of  the  pipe¬ 
line  processed  data,  which  is  taken  to  be  <20%  for  SL  and  LL 
(Infrared  Spectrograph  Data  Handbook).  This  level  of  possible 
flux  uncertainty  will  not  affect  our  main  discussions  significantly. 
In  fact,  in  the  majority  of  ULIRGs  with  IRAS  detection  at  25  pm, 
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Fig.  2. — Schematic  diagram  of  a  strong  dust  temperature  gradient  around  a 
centrally  concentrated  energy  source  ( central  small  filled  circle).  Inner  dust 
close  to  the  central  energy  source  has  higher  temperature  than  dust  far  away 
from  the  energy  source.  Assuming  approximately  blackbody  radiation,  dust  at 
1000  K  dust,  close  to  the  innermost  dust  sublimation  radius,  shows  an  emission 
peak  at  2rest  ~  3  /im.  More  distant  300  K  dust  and  even  150  K  dust  have  emis¬ 
sion  peaks  at  2rest  ~  10  and  ~20  fin i,  respectively.  This  is  a  very  simplified  pic¬ 
ture,  and  in  an  actual  case,  dust  with  other  temperatures  can  also  contribute  to  the 
observed  flux  at  each  wavelength.  For  any  reasonable  dust  radial  density  distribu¬ 
tion,  it  is  always  true  that  the  contribution  from  inner,  hotter  dust  emission  is 
larger  at  a  shorter  wavelength  in  this  wavelength  range.  When  we  look  at  the 
energy  source  from  the  right  side,  an  observed  flux  is  dominated  by  dust  at  the 
foreground  side  {right),  rather  than  dust  with  the  same  temperature  at  the  distant 
side  {left),  due  to  flux  attenuation,  particularly  for  outer  cooler  dust.  Thus,  the 
relation  Ay(3  Mm)  >  Ay{10  fin l)  >  Ay{ 20  fim)  should  hold.  Dust  extinction  to¬ 
ward  the  3  fim  continuum-emitting  region,  Ay {3  fim),  can  be  estimated  from  the 
optical  depths  of  absorption  features  at  3.4  fin i  due  to  bare  carbonaceous  dust  and 
at  3.1  fin i  due  to  ice-covered  dust  (Imanishi  et  al.  2001,  2006a;  Imanishi  2006). 
The  extinction  toward  the  10  and  20  fim  continuum-emitting  regions,  Ay {10  fim) 
and  Ay  {20  fim),  can  be  estimated  from  the  optical  depths  of  9.7  and  18  fim  silicate 
dust  absorption  features  in  Spitzer  IRS  spectra,  respectively.  Although  a  spherical 
geometry  is  illustrated  here  for  simplicity,  the  essence  is  unchanged  as  long  as  a  cen¬ 
trally  concentrated  compact  energy  source  is  obscured  by  dust  along  our  line  of  sight. 


the  Spitzer  IRS  25  pm  flux  agrees  within  30%  to  the  IRAS  25  pm 
data,  and  the  Spitzer  IRS  flux  is  generally  smaller.  In  a  small  frac¬ 
tion  of  ULIRGs,  Spitzer  IRS  25  pm  flux  is  =50%  of  IRAS  25  pm 
data,  possibly  because  extended  dim  emission  is  important.  For 
ULIRGs  with  IRAS  nondetection  at  25  pm,  the  measured  Spizter 
IRS  flux  at  25  pm  is  always  smaller  than  the  upper  limit  of  the 
IRAS  data. 

For  a  fraction  of  the  ULIRGs,  slight  flux  discrepancies  between 
SL1  and  LL2  were  discernible,  ranging  from  10%  to  40%.  When 
the  discrepancy  is  seen,  the  SL1  flux  (3.7"  wide  slit)  is  usually 
smaller  than  the  LL2  flux  (10.5"),  suggesting  the  presence  of  ex¬ 
tended  emission  in  these  ULIRGs.  In  these  cases  we  adjusted  the 
smaller  SL 1  (and  SL2)  flux  to  match  the  larger  LL2  flux. 

The  main  part  of  the  9.7  pm  silicate  absorption  feature  extends 
Arest  =  8-13  pm  in  the  rest  wavelength  (e.g.,  Chiar  &  Tielens 
2006).  For  ULIRGs  at  z  >  0. 1 1,  the  feature  is  not  fully  covered 
with  SL1  spectra  at  A0bs  =  7.4-14.5  pm.  Even  for  ULIRGs  at 
z  <  0 . 1 1 ,  the  continuum  data  at  lrest  =  13-14  ^m,  which  are  re¬ 
quired  to  estimate  the  strength  of  the  9.7  pm  silicate  absorption, 
fall  beyond  Aobs  >  14.5  pm,  outside  the  wavelength  range  of  SL  I . 
Thus,  for  the  majority  of  the  observed  ULIRGs,  LL2  spectra  are 
essential  to  measure  the  optical  depth  of  the  9.7  pm  silicate  ab¬ 
sorption  feature  with  a  small  ambiguity. 

Appropriate  spectral  binning  with  2  or  4  pixels  was  applied  to 
reduce  the  scatter  of  data  points  at  SL2  (5. 2-7. 7  pm)  for  some 
faint  ULIRGs,  and  at  A0bs  ~  10  pm  for  ULIRGs  displaying  very 
strong  9.7  pm  silicate  absorption  features. 


4.  RESULT 

Figures  3  and  4  present  low-resolution  5-35  pm  spectra  of 
ULIRGs  classified  optically  as  LINERs  and  H  n  regions,  respec¬ 
tively.  For  most  of  the  sources,  full  5-35  pm  spectra  are  shown 
here  for  the  first  time.  Although  spectra  of  the  five  ULIRGs  in 
the  GTO  program  (Arp  220,  IRAS  08572+3915,  12112+0305, 
14348— 1447,  and  2249 1  —  1 808)  were  presented  elsewhere  (Spoon 
et  al.  2006;  Annus  et  al.  2007),  we  include  their  spectra  to  com¬ 
prise  a  complete  sample,  to  which  our  detailed  and  uniform  anal¬ 
ysis  is  applied.  For  these  ULIRGs,  their  and  our  spectra  are,  in 
overall,  consistent  with  each  other. 

The  spectra  in  Figures  3  and  4  are  suitable  for  displaying  the 
properties  of  the  9.7  and  1 8  pm  silicate  dust  absorption  features. 
However,  they  are  not  very  useful  for  PAH  emission  features. 
Figure  5  presents  zoom-in  spectra  at  l0bS  =  5.2-14.5  pm  of  op¬ 
tical  LINER  and  H  n  region  ULIRGs  to  better  examine  the  prop¬ 
erties  of  the  PAH  emission  features. 

4. 1 .  PAH  Emission 

The  majority  of  ULIRGs  in  Figure  5  show  clearly  detectable 
PAH  emission  features  at  2rest  =  6.2,  7.7,  and  11.3  pm.  A  frac¬ 
tion  of  the  ULIRGs  also  display  the  PAH  emission  feature  at 
^rest=  8.6  pm.  Emission  at  Arest  =  12.8  pm  is  also  detected 
(Figs.  3  and  4),  but  it  can  be  attributed  to  both  the  12.8  pm  PAH 
and  [Ne  ii]  emission  lines,  which  are  not  clearly  resolved  in  our 
low-resolution  spectra.  We  thus  will  not  discuss  the  12.8  pm 
PAH  emission  feature  further. 

The  6.2  pm  PAH  emission  feature  is  attributed  to  an  aromatic 
C-C  stretching  fundamental  (Allamandola  et  al.  1 989).  It  is  rela¬ 
tively  weak  but  is  isolated  both  from  other  PAH  emission  and 
from  the  strong  9.7  pm  silicate  dust  absorption  feature.  Although 
absorption  at  6.0  pm  due  to  H20  ice  and  at  6.85  pm  by  hydro¬ 
genated  amorphous  carbon  (HAC)  are  present  (Spoon  et  al.  2002), 
their  relatively  weak  nature,  compared  to  the  9.7  pm  silicate  ab¬ 
sorption,  makes  it  possible  to  measure  the  6.2  pm  PAH  emission 
strength  with  little  ambiguity.  For  this  reason,  the  6.2  pm  feature 
has  often  been  used  to  measure  the  PAH  emission  strength  in  a 
UL1RG  (  Fischer  2000;  Spoon  et  al.  2007).  To  estimate  the  strength 
of  the  6.2  pm  PAH  emission  feature,  we  adopt  a  linear  continuum 
determined  from  data  at  Arest  =  6.1  and  6.45  pm,  which  is  shown 
as  solid  lines  in  Figure  5. 

The  11.3  pm  PAH  emission  is  due  to  an  out-of-plane  aromatic 
C-H  bending  vibration  (Allamandola  et  al.  1989).  Since  it  is 
situated  close  to  the  peak  of  the  strong  9.7  pm  silicate  absorption, 
its  flux  attenuation  can  be  severe  if  the  PAH-emitting  regions  are 
highly  dust  obscured.  Thus,  its  flux  relative  to  other  PAH  fea¬ 
tures  can  be  a  good  indicator  of  the  dust  obscuration  of  the  PAH- 
emitting  regions.  To  estimate  the  1 1 .3  pm  PAH  strength,  we  adopt 
a  linear  continuum  connecting  data  points  at  Arest  =11.0  and 
11.6  pm,  which  is  shown  as  solid  lines  in  Figure  5. 

The  7.7  pm  PAH  feature  is  due  to  aromatic  C-C  stretching 
vibrations  (Allamandola  et  al.  1989).  Although  it  is  the  strongest 
PAH  emission  feature,  its  substantial  wavelength  overlap  with 
the  broad,  strong  9.7  pm  silicate  dust  absorption  feature  makes  it 
very  difficult  to  distinguish  between  7.7  pm  PAH  emission  and 

9.7  pm  silicate  absorption,  particularly  for  sources  showing  strong 
silicate  absorption,  such  as  the  majority  of  ULIRGs.  To  reduce  the 
effects  of  the  9.7  pm  absorption  feature,  we  adopt  an  excess  above 
a  linear  continuum  determined  from  data  at  Arest  =  7.3  and  8.1  pm 
(shown  as  solid  lines  in  Fig.  5)  as  the  7.7  pm  PAH  emission.  Since 
many  previous  papers  (e.g.,  Genzel  et  al.  1998)  measured  the 

7.7  pm  PAH  fluxes  based  on  different  definitions,  readers  must 
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TABLE  1 

Observed  ULIRGs  at  z  <  0.15  and  Their  IRAS-b ased  Infrared  Emission  Properties 


Object 

(1) 

Redshift 

(2) 

f\2 

(Jy) 

(3) 

fl5 

(Jy) 

(4) 

k  0 
(Jy) 

(5) 

/ioo 

(Jy) 

(6) 

log  Hr 
(i0) 

(7) 

fislfeo 

(8) 

Optical 

Class 

(9) 

IRAS  00188-0856 . 

0.128 

<0.12 

0.37 

2.59 

3.40 

12.3 

0.14  (C) 

LINER 

IRAS  00482-2721 . 

0.129 

<0.10 

<0.18 

1.13 

1.84 

12.0 

<0.16  (C) 

LINER 

IRAS  03250+1606 . 

0.129 

<0.10 

<0.15 

1.38 

1.77 

12.1 

<0.11  (C) 

LINER 

IRAS  04103-2838 . 

0.118 

0.08 

0.54 

1.82 

1.71 

12.2 

0.30  (W) 

LINER 

IRAS  08572+3915 . 

0.058 

0.32 

1.70 

7.43 

4.59 

12.1 

0.23  (W) 

LINER 

IRAS  09039+0503 . 

0.125 

0.07 

0.12 

1.48 

2.06 

12.1 

0.08  (C) 

LINER 

IRAS  09116+0334 . 

0.146 

<0.09 

<0.14 

1.09 

1.82 

12.1 

<0.13  (C) 

LINER 

IRAS  09539+0857 . 

0.129 

<0.15 

<0.15 

1.44 

1.04 

12.0 

<0.11  (C) 

LINER 

IRAS  10378+1108 . 

0.136 

<0.11 

0.24 

2.28 

1.82 

12.3 

0.11  (C) 

LINER 

IRAS  10485-1447 . 

0.133 

<0.11 

0.25 

1.73 

1.66 

12.2 

0.14  (C) 

LINER 

IRAS  10494+4424 . 

0.092 

<0.12 

0.16 

3.53 

5.41 

12.1 

0.05  (C) 

LINER 

IRAS  11095-0238 . 

0.106 

0.06 

0.42 

3.25 

2.53 

12.2 

0.13  (C) 

LINER 

IRAS  11130-2659 . 

0.136 

<0.09 

0.20 

1.21 

1.24 

12.1 

0.17  (C) 

LINER 

IRAS  12112+0305 . 

0.073 

0.12 

0.51 

8.50 

9.98 

12.3 

0.06  (C) 

LINER 

IRAS  12127-1412 . 

0.133 

<0.13 

0.24 

1.54 

1.13 

12.1 

0.16  (C) 

LINER 

IRAS  12359-0725 . 

0.138 

0.09 

0.15 

1.33 

1.12 

12.1 

0.11  (C) 

LINER 

IRAS  13335-2612 . 

0.125 

<0.13 

<0.14 

1.40 

2.10 

12.1 

<0.10  (C) 

LINER 

IRAS  14252-1550 . 

0.149 

<0.09 

<0.23 

1.15 

1.86 

12.2 

<0.20  (C) 

LINER 

IRAS  14348-1447 . 

0.083 

0.07 

0.49 

6.87 

7.07 

12.3 

0.07  (C) 

LINER 

IRAS  15327+2340  (Arp  220) . 

0.018 

0.48 

7.92 

103.33 

112.40 

12.1 

0.08  (C) 

LINER 

IRAS  16090-0139 . 

0.134 

0.09 

0.26 

3.61 

4.87 

12.5 

0.07  (C) 

LINER 

IRAS  16468+5200 . 

0.150 

<0.06 

0.10 

1.01 

1.04 

12.1 

0.10  (C) 

LINER 

IRAS  16487+5447 . 

0.104 

<0.07 

0.20 

2.88 

3.07 

12.1 

0.07  (C) 

LINER 

IRAS  17028+5817 . 

0.106 

<0.06 

0.10 

2.43 

3.91 

12.1 

0.04  (C) 

LINER 

IRAS  17044+6720 . 

0.135 

<0.07 

0.36 

1.28 

0.98 

12.1 

0.28  (W) 

LINER 

IRAS  21329-2346 . 

0.125 

0.05 

0.12 

1.65 

2.22 

12.1 

0.07  (C) 

LINER 

IRAS  23234+0946 . 

0.128 

<0.06 

0.08 

1.56 

2.11 

12.1 

0.05  (C) 

LINER 

IRAS  23327+2913 . 

0.107 

<0.06 

0.22 

2.10 

2.81 

12.1 

0.10  (C) 

LINER 

IRAS  00091-0738 . 

0.118 

<0.07 

0.22 

2.63 

2.52 

12.2 

0.08  (C) 

H  ii 

IRAS  00456-2904 . 

0.110 

<0.08 

0.14 

2.60 

3.38 

12.2 

0.05  (C) 

H  ii 

IRAS  01004-2237 . 

0.118 

0.11 

0.66 

2.29 

1.79 

12.3 

0.29  (W) 

H  ii 

IRAS  01166-0844 . 

0.118 

0.07 

0.17 

1.74 

1.42 

12.1 

0.10  (C) 

H  ii 

IRAS  01298-0744 . 

0.136 

<0.12 

0.19 

2.47 

2.08 

12.3 

0.08  (C) 

H  ii 

IRAS  01569-2939 . 

0.141 

<0.11 

0.14 

1.73 

1.51 

12.2 

0.08  (C) 

H  ii 

IRAS  02411+0353 . 

0.144 

<0.08 

0.22 

1.37 

1.95 

12.2 

0.16  (C) 

H  ii 

IRAS  10190+1322 . 

0.077 

<0.07 

0.38 

3.33 

5.57 

12.0 

0.11  (C) 

H  ii 

IRAS  11387+4116 . 

0.149 

0.12 

<0.14 

1.02 

1.51 

12.2 

<0.14  (C) 

H  ii 

IRAS  11506+1331  . 

0.127 

<0.10 

0.19 

2.58 

3.32 

12.3 

0.07  (C) 

H  ii 

IRAS  13509+0442 . 

0.136 

0.10 

<0.23 

1.56 

2.53 

12.2 

<0.15  (C) 

H  ii 

IRAS  13539+2920 . 

0.108 

<0.09 

0.12 

1.83 

2.73 

12.0 

0.07  (C) 

H  ii 

IRAS  14060+2919 . 

0.117 

<0.10 

0.14 

1.61 

2.42 

12.1 

0.09  (C) 

H  ii 

IRAS  15206+3342 . 

0.125 

0.08 

0.35 

1.77 

1.89 

12.2 

0.20  (C) 

H  ii 

IRAS  15225+2350 . 

0.139 

<0.07 

0.18 

1.30 

1.48 

12.1 

0.14(C) 

H  ii 

IRAS  16474+3430 . 

0.111 

<0.13 

0.20 

2.27 

2.88 

12.1 

0.09  (C) 

H  ii 

IRAS  20414-1651 . 

0.086 

<0.65 

0.35 

4.36 

5.25 

12.2 

0.08  (C) 

H  ii 

IRAS  21208-0519 . 

0.130 

<0.09 

<0.15 

1.17 

1.66 

12.0 

<0.13  (C) 

H  ii 

IRAS  22206-2715 . 

0.132 

<0.10 

<0.16 

1.75 

2.33 

12.1 

<0.09  (C) 

H  ii 

IRAS  22491-1808 . 

0.076 

0.05 

0.55 

5.44 

4.45 

12.1 

0.10  (C) 

H  ii 

Notes. — Col.  (1):  Object  name.  Col.  (2):  Redshift.  Cols.  (3) — (6): /12,  and  /100  are  IRAS  fluxes  at  12,25, 60,  and  100  /im,  respectively,  taken  from  Kim  &  Sanders 

(1998).  Col.  (7):  Decimal  logarithm  of  infrared  (8-1000  ^m)  luminosity  in  units  of  solar  luminosity  ( LQ ),  calculated  with  =  2.1  x  1039  x  Z)(Mpc)2  x  (13.48  x 
fi2  +  5.16  x  fz 5  +  2.58  x  0  +  /loo)  ergs  s-1  (Sanders  &  Mirabel  1996).  Since  the  calculation  is  based  on  our  adopted  cosmology,  the  infrared  luminosities  slightly  (<  10%) 
differ  from  the  values  shown  in  Kim  &  Sanders  (1998,  their  Table  1,  col.  [15]).  For  sources  that  have  upper  limits  in  some  IRAS  bands,  we  can  derive  upper  and  lower  limits 
of  the  infrared  luminosity  by  assuming  that  the  actual  flux  is  the  IRAS  upper  limit  and  zero  value,  respectively.  The  difference  of  the  upper  and  lower  values  is  usually  very 
small,  less  than  0.2  dex.  We  assume  that  the  infrared  luminosity  is  the  average  of  these  values.  Col.  (8):  IRAS  25-60  fim  flux  ratio.  ULIRGs  with  fis/feo  <  0-2  and  >0.2 
are  classified  as  cool  and  warm  sources  (denoted  as  “C”  and  “W”),  respectively  (Sanders  et  al.  1988b).  Col.  (9):  Optical  spectral  classification  by  Veilleux  et  al.  (1999a). 


be  careful  when  our  value  is  compared  with  other  estimates  in 
the  literature.  In  this  paper,  although  the  7.7  /im  PAH  strength  is 
shown  for  reference,  it  will  not  be  used  when  possible  uncertainty 
coming  from  the  effects  of  the  strong  9.7  /im  silicate  feature  may 
complicate  quantitative  discussions. 


The  8.6  /im  PAH  emission  due  to  in-plane  aromatic  C-H  bend¬ 
ing  vibration  (Allamandola  et  al.  1989)  is  not  used  in  this  paper, 
because  it  is  weak  and  overlaps  with  the  9.7  /im  silicate  absorp¬ 
tion,  which  make  its  strength  highly  uncertain.  Other  weak  PAH 
emission  features  may  be  present  (e.g.,  Arest  =  5.7  /im;  Spoon 


TABLE  2 

Spitzer  IRS  Observing  Log 

Integration  Time  (s) 


Object  PID  Date  (UT)  SL2  SL1  LL2  LL1 

(1)  (2)  (3)  (4)  (5)  (6)  (7) 


IRAS  00188-0856 .  105  2003  Dec  17  240  240  180  180 

IRAS  00482-2721  .  3187  +  2306  2005  Jul  7  +  2004  Jul  1 8  240  240  240  240 

IRAS  03250+1606 .  3187  +  2306  2005  Feb  11  240  240  240  240 

IRAS  04103-2838 .  3187  +  2306  2005  Feb  10  +  2004  Aug  10  240  240  240  240 

IRAS  08572+3915 .  105  2004  Apr  15  84  84  84  84 

IRAS  09039+0503 .  3187  +  2306  2005  Apr  1 8  +  2004  Nov  1 7  240  240  240  240 

IRAS  09116+0334 .  2306  2005  Apr  21  240  240  240  240 

IRAS  09539+0857 .  3187  +  2306  2005  Jun  5  +  2005  May  3 1  240  240  240  240 

IRAS  10378+1108 .  105  2005  Jun  8  240  240  180  180 

IRAS  10485-1447 .  3187  +  2306  2005  May  23  +  2005  Jan  4  240  240  480a  480a 

IRAS  10494+4424 .  2306  2004  Nov  17  240  240  240  240 

IRAS  11095-0238 .  105  2005  Jun  7  240  120  240  120 

IRAS  11130-2659 .  2306  2005  Jul  12  240  240  240  240 

IRAS  12112+0305 .  105  2004  Jan  4  84  84  120  120 

IRAS  12127-1412 .  3187  +  2306  2005  Jun  30  240  240  240  240 

IRAS  12359-0725 .  2306  2005  Jun  30  240  240  240  240 

IRAS  13335-2612 .  3187  +  2306  2005  Feb  15  +  2004  Jul  17  240  240  240  240 

IRAS  14252-1550 .  2306  2004  Jul  17  240  240  240  240 

IRAS  14348-1447 .  105  2004  Feb  7  120  120  120  120 

Arp  220 .  105  2004  Feb  29  84  84  60  60 

IRAS  16090-0139 .  105  2005  Aug  6  120  120  180  180 

IRAS  16468+5200 .  2306  2004  Jul  14  480b  480b  480a  480a 

IRAS  16487+5447 .  2306  2004  Jul  17  240  240  240  240 

IRAS  17028+5817 .  2306  2004  Jul  17  240  240  240  240 

IRAS  17044+6720 .  2306  2004  Jul  17  240  240  240  240 

IRAS  21329-2346 .  3187  +  2306  2004  Nov  16  +  2004  Oct  23  240  240  240  240 

IRAS  23234+0946 .  3187  +  2306  2004  Dec  13  +  2004  Dec  8  240  240  240  240 

IRAS  23327+2913 .  2306  2004  Dec  8  240  240  240  240 

IRAS  00091-0738 .  3187  +  2306  2005  Jun  30  +  2004  Dec  8  240  240  240  240 

IRAS  00456-2904 .  3187  +  2306  2005  Jul  14  +  2004  Jul  17  240  240  240  240 

IRAS  01004-2237 .  105  2004  Jan  4  120  120  120  120 

IRAS  01166-0844 .  3187  +  2306  2005  Jan  3  +  2005  Jan  15  +  2004  Jul  17  480c  480°  480a  480a 

IRAS  01298-0744 .  105  2005  Jul  14  240  240  120  120 

IRAS  01569-2939 .  2306  2004  Jul  18  240  240  240  240 

IRAS  02411+0353 .  2306  2005  Jan  14  240  240  240  240 

IRAS  10190+1322 .  3187  +  2306  2005  May  22  +  2004  Dec  12  240  +  240d  240  +  240d  480a  480a 

IRAS  11387+4116 .  2306  2005  Jan  11  240  240  240  240 

IRAS  11506+1331  .  3187  +  2306  2005  May  25  +  2005  Jan  3  240  240  240  240 

IRAS  13509+0442 .  2306  2004  Jul  17  240  240  240  240 

IRAS  13539+2920 .  2306  2005  Feb  7  240  240  240  240 

IRAS  14060+2919 .  2306  2004  Jul  16  240  240  240  240 

IRAS  15206+3342 .  105  2004  Jun  24  120  120  180  180 

IRAS  15225+2350 .  2306  2005  Feb  7  240  240  240  240 

IRAS  16474+3430 .  2306  2004  Jul  14  240  240  240  240 

IRAS  20414-1651  .  105  2004  May  14  120  120  120  120 

IRAS  21208-0519 .  3187  +  2306  2004  Nov  13  +  2004  Oct  22  240e  240e  480a  480a 

IRAS  22206-2715 .  3187  +  2306  2004  Nov  15  +  2004  Oct  25  240  240  480a  480a 

IRAS  22491-1808 .  105  2004  Jun  24  120  120  120  120 


Notes. — Col.  (1):  Object  name.  Col.  (2):  PID  number:  105  (PI:  J.  Houck),  2306  (PI:  M.  Imanishi),  and  3187  (PI:  S.  Veilleux).  For  sources  denoted  with  “3187  + 
2306,”  SL  and  LL  spectra  were  taken  by  observations  with  PIDs  =  3187  and  2306,  respectively.  Col.  (3):  Observing  date  in  UT.  Col.  (4):  Net  on-source  integration  time 
for  SL2  spectroscopy  in  seconds.  Col.  (5):  Net  on-source  integration  time  for  SL1  spectroscopy  in  seconds.  Col.  (6):  Net  on-source  integration  time  for  LL2  spec¬ 
troscopy  in  seconds.  Col.  (7):  Net  on-source  integration  time  for  LL1  spectroscopy  in  seconds. 

a  Although  separate  spectroscopy  of  double  nuclei  (Kim  et  al.  2002)  was  proposed,  the  same  coordinate  was  observed  twice,  probably  because  IRS  blue  peak  up 
camera  (13-18.5  pm)  pointed  to  the  same  position. 

b  We  combined  SL  spectra  of  double  nuclei  with  ~3"  separation  (Kim  et  al.  2002),  because  both  nuclei  are  very  faint.  The  combined  SL1  spectrum  is  smoothly 
connected  to  the  LL2  spectrum. 

c  We  combined  SL  spectra  of  double  nuclei  with  ~5"  separation  (Kim  et  al.  2002),  because  the  fainter  northern  nucleus  is  very  faint.  The  combined  SL1  spectrum  is 
smoothly  connected  to  the  LL2  spectrum. 

d  SL  spectra  of  the  eastern  and  western  nuclei  with  ~4"  separation  (Kim  et  al.  2002)  are  analyzed  separately,  because  both  nuclei  have  similar  flux  levels.  The 
combined  SL1  spectrum  is  smoothly  connected  to  the  LL2  spectrum. 

e  Only  the  brighter  northern  nucleus  (Kim  et  al.  2002)  was  analyzed,  because  (1)  the  fainter  southern  nucleus  with  ~7"  separation  is  very  faint,  (2)  it  may  not  be  fully 
covered  with  the  LL  spectrum,  and  (3)  the  SL1  spectrum  of  the  northern  nucleus  is  smoothly  connected  to  the  LL2  spectrum. 
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Fig.  3. — Spectra  of  ULIRGs  classified  optically  as  LINERs.  The  abscissa  and  ordinate  are,  respectively,  the  observed  wavelength  in  /im  and  flux  Fv  in  Jy,  both  shown 
as  decimal  logarithmic  scale.  For  all  objects,  the  ratio  of  the  uppermost  to  lowermost  scale  in  the  ordinate  is  a  factor  of  1000,  to  show  the  variation  of  the  overall  spectral 
energy  distribution.  Dashed  line :  power-law  continuum  determined  from  data  points  at  /.rest  =  5.6  and  7.1  fii n,  and  20bS  =  34—35  /im.  Dotted  line',  power-law  con¬ 
tinuum  determined  from  data  points  at  yLrest  =  7.1  and  14.2  fim  for  9.7  /im  silicate  absorption,  and  at  2rest  =  14.2  and  24  /im  for  18  /im  silicate  absorption.  Dash-dotted 
line :  power-law  continuum  determined  from  data  points  at  /.rest  =  14.2  and  29  /im,  used  only  for  ULIRGs  showing  strong  18  /im  silicate  dust  absorption  features. 


et  al.  2002)  but  are  not  discussed  in  detail  because  no  new  mean¬ 
ingful  information  is  expected  to  be  obtained. 

In  summary,  the  6.2  and  11.3  pm  PAH  emission  features  are 
the  main  focus  of  our  discussion,  and  the  strength  of  the  7.7  /im 
PAH  feature  based  on  our  own  definition  is  shown  just  for  refer¬ 
ence.  We  fit  the  PAH  emission  features  with  Gaussian  profiles, 
which  reproduce  the  observed  data  reasonably  well.  The  observed 
rest- frame  equivalent  widths  and  luminosities  of  the  6.2,  7.7,  and 
11.3  //m  PAH  emission  features,  based  on  our  adopted  continuum 
levels,  are  summarized  in  Table  3.  The  uncertainties  coming  from 
the  fittings  are  unlikely  to  exceed  30%,  so  once  the  continuum 
levels  are  fixed,  the  total  uncertainties  are  of  this  order. 

4.2.  Silicate  Absorption 

Based  on  our  adopted  continuum  levels,  the  9.7  pm  silicate 
absorption  feature  is  found  in  almost  all  ULIRGs,  and  1 8  pm 
silicate  absorption  is  also  seen,  particularly  in  ULIRGs  with 
strong  9.7  pm  absorption  (Figs.  3  and  4).  Silicate  dust  grains  show 
two  strong  absorption  peaks  centered  at  2rest  =  9.7  and  18.5  pm. 
Outside  these  peaks,  a  weak  absorption  tail  extends  from  2rest  =  8 


to  >30  pm  (Chiar  &  Helens  2006;  see  also  Fig.  6 a  in  this  paper). 
Thus,  all  wavelengths  between  Arest  =  8  and  30  pm  show  excess 
absorption,  compared  to  the  extrapolation  of  dust  extinction  curve 
from  wavelengths  outside  Arest  =  8-30  pm.  Hence,  to  accurately 
estimate  the  optical  depths  of  the  9.7  pm  (r9  7)  and  18  pm  (r18) 
silicate  absorption  features,  a  continuum  level  must  be  determined 
using  data  at  Arest  <  8  and  >30  pm.  At  /lrest  =  6-8  pm.  PAH  emis¬ 
sion  can  contaminate  the  continuum  level,  so  that  in  practice  data 
points  at  Arest  ~  5.6  pm  are  used  at  the  shorter  wavelength  side. 
Data  points  at  lrest  ~  7. 1  pm  are  also  less  affected  by  PAH  emis¬ 
sion  and  absorption  features  (Spoon  et  al.  2006)  and  so  can  be 
used  as  additional  continuum  points.  Table  4  summarizes  the 
optical  depths  of  the  9.7  pm  silicate  absorption  (t9  7),  against 
power-law  continuum  levels  determined  from  data  points  at 
lrest  =  5.6  and  7.1  pm,  and  l0bs  =  34-35  pm  (the  longest  part 
of  the  IRS  spectra)  shown  as  the  dashed  lines  in  Figures  3  and 
4.  A  power-law  continuum  is  assumed  here,  because  the  intrinsic 
unabsorbed  5-35  /mi  spectrum  of  AGN-heated  hot  dust  emis¬ 
sion  is  roughly  approximated  by  a  power-law  form  (Weedman 
et  al.  2005). 
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Fig.  3  —  Continued 


The  wavelengths  of  these  data  points  used  for  the  continuum 
determination  are,  however,  far  away  from  the  9.7  and  18  pm 
silicate  peaks.  The  optical  depth  values  of  the  9.7  and  1 8  pm  sili¬ 
cate  absoiption  features  can  vary  depending  on  the  assumptions 
of  the  continuum  data  points  as  well  as  the  continuum  shape 
(power-law,  linear,  spline,  or  polynomial,  etc.).  In  addition,  the 
power-law  continua  do  not  necessarily  provide  perfect  fits  to 
observed  data  points  outside  the  silicate  features.  In  the  spec¬ 
trum  of  the  weakly  obscured  ULIRG1RAS  01004—2237  (Fig.  4, 
third  plot),  data  points  at  A0bs  =  13-18  pm  are  clearly  above  the 
power-law  continuum  ( dashed  line),  indicating  that  the  power- 
law  continuum  only  roughly  represents  the  overall  continuum. 
The  relative  ambiguity  of  the  T97  value  is  small  because  the  9.7  pm 
silicate  feature  is  intrinsically  strong.  However,  the  optical  depth 
of  the  much  weaker  1 8  pm  silicate  absorption  feature  can  be 
highly  uncertain.  Consequently,  discussions  of  the  presence  or 
absence  of  strong  dust  temperature  gradients,  based  on  the  rel¬ 
ative  strength  of  the  9.7  and  18  pm  silicate  dust  absorption  fea¬ 
tures  (Fig.  2),  could  also  be  highly  uncertain,  if  we  use  T9.7  and 
rig  as  defined  above. 

For  this  reason,  we  define  Tg  7  and  t[ 8 .  The  r9'  7  value  is  the  op¬ 
tical  depth  of  the  9.7  pm  silicate  feature  against  a  power-law  con¬ 


tinuum  determined  from  data  points  at  Arest  =  7.1  and  14.2  pm. 
The  ri's  value  is  the  optical  depth  of  the  1 8  pm  silicate  feature 
against  a  power-law  continuum  determined  from  data  points  at 
Arest  =  14.2  and  24  pm.  These  continua  are  shown  as  dotted  lines 
in  Figures  3  and  4.  With  these  new  definitions,  the  ambiguities 
are  reduced,  because  the  continuum  levels  are  determined  using 
data  points  just  outside  the  9.7  and  1 8  pm  features,  closer  to  the 
absorption  peaks.  The  rj  7  values  are  shown  for  all  ULIRGs  in 
Table  4  (col.  [3]).  The  Tls  values  are  also  shown  in  Table  4  (col.  [4]) 
for  ULIRGs,  where  the  1 8  pm  feature  is  clearly  seen  in  absorption. 

The  Tls  values  are  appropriate  to  estimate  the  strength  of  the 
1 8  pm  silicate  feature,  if  the  1 8  pm  profile  is  similar  to  that  ob¬ 
served  in  the  ULIRG  IRAS  08572+3915  (Spoon  et  al.  2006). 
However,  the  18  pm  silicate  profile  of  the  Galactic,  highly  ob¬ 
scured  star  GCS3-I  extends  far  beyond  Arest  ~  24  pm  and  reaches 
to  Arest  ~  29  pm  (Chiar  &  Helens  2006;  see  also  Fig.  6 a  in  this 
paper).  Thus,  we  also  define  Tj'g,  which  is  the  optical  depth  of  the 
18  pm  feature  against  a  power-law  continuum  determined  from 
data  points  at  Arest  =  14.2  and  29  pm  (shown  as  dash-dotted  lines 
in  Figs.  3  and  4).  The  r[g  values  are  also  summarized  in  Table  4 
(col.  [5])  for  ULIRGs  with  clearly  detectable  18  pm  silicate 
absorption. 


Fig.  3  —  Continued 


Fig.  4. — Spectra  of  ULIRGs  classified  optically  as  H  n  regions,  shown  in  the  same  way  as  Fig.  3.  Symbols  are  also  the  same  as  Fig.  3.  For  IRAS  10190+1322,  SL 
spectra  of  the  eastern  and  western  nuclei  are  combined. 
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Fig.  4 — Continued 


4.3.  Ice  Absorption 

A  significant  fraction  of  ULIRGs  display  dips  at  the  shorter 
wavelength  part  of  the  6.2  pm  PAH  emission  features.  To  better 
visualize  the  dips,  Figure  7  presents  zoom-in  spectra  at  A0bs  = 
5.2-9  pm  for  ULIRGs  showing  clear  or  possible  dips.  Spectra  of 
two  ULIRGs  without  discernible  dips  (IRAS  17044+6720  and 
00456—2904)  are  also  shown  as  examples  of  no  clear  detection. 
We  ascribe  the  dips  to  the  6.0  pm  H20  absorption  feature  (bend¬ 
ing  mode),  as  has  been  found  in  both  Galactic  obscured  stars 
(Gibb  et  al.  2000;  Keane  et  al.  2001)  and  obscured  ULIRGs  (Spoon 
et  al.  2002).  Since  the  6.0  pm  ice  absorption  feature  is  sand¬ 
wiched  by  PAH  emission  features  at  Arest  =  5.7  and  6.2  pm,  it  is 
sometimes  difficult  to  distinguish  whether  the  dip  is  due  to  ice 
absorption  or  to  PAH  emission  features.  For  ULIRGs  with  clearly 
detectable  6.0  pm  H20  ice  absorption  features,  their  observed 
optical  depths  (r6  0)  are  summarized  in  Table  5. 

For  ULIRGs  showing  clear  ice  (and  HAC)  absorption  features, 
the  continuum  flux  levels  just  outside  of  the  6.2  pm  PAH  emission 
are  reduced.  For  these  ULIRGs,  approximate  EW6  2PAH  values 
after  these  corrections  are  also  shown  in  Table  3 . 

5.  DISCUSSION 

5.1.  Magnitudes  of  Detected  Starbursts 

PAH  emission  luminosities  can  roughly  probe  the  absolute 
magnitudes  of  PAH-emitting  normal  starbursts,  if  the  extinction 
of  the  starbursts  is  say  less  than  about  20-30  mag  in  Ay,  because 
this  amount  of  Ay  produces  little  flux  attenuation  (<1  mag)  at 
2.rest  >  5  pm,  aside  from  the  strong  9.7  pm  silicate  absorption 
feature  (Rieke  &  Lebofsky  1985;  Lutz  et  al.  1996). 

Table  3  (cols.  [9]  and  [10])  tabulates  the  6.2  pm  PAH-to-infrared 
luminosity  ratio,  L6  2 pah/^ir,  and  the  1 1 .3  pm  PAH-to-infrared  lumi¬ 
nosity  ratio,  Zn  3PAH/LIR.  The  ratios  in  normal  starburst  galaxies  with 


modest  dust  obscuration  (A  v  <  20-30  mag)  are  estimated  to  be 
7-6.2pah/7.ir  ~  3.4  x  10~3  (Peetersetal.  2004)  and Lujpah/Yir  ~ 
1 .4  x  1 0-3  (Soifer  et  al.  2002).  Figures  8 a  and  8 b  display  the  PAH- 
to-infrared  luminosity  ratio  ( ordinate )  and  PAH  equivalent  width 
( abscissa )  for  the  6.2  and  11.3  pm  PAH  features,  respectively.  The 
observed  T6.2pah/7-ir  and  Lh.3Pah/7.ir  ratios  in  the  majority  of 
the  observed  ULIRGs  are  systematically  lower  than  those  seen  in 
modestly  obscured  normal  starburst  galaxies,  suggesting  that 
the  detected  starburst  activity  as  measured  by  the  PAH  emission 
in  ULIRGs  is  energetically  insignificant.  The  dominant  energy 
sources  of  ULIRGs  must  therefore  be  either  (1)  highly  obscured 
(Ay  »  30  mag)  starbursts,  where  the  emitted  PAH  luminosities 
are  substantially  attenuated  by  dust  extinction,  or  (2)  buried  AGNs, 
which  produce  intrinsically  weak  PAH  emission. 

5.2.  Candidates  of  ULIRGs  Possessing  Powerful  Buried  AGNs 

5.2.1.  Low  Equivalent  Widths  of  PAH  Emission  Features 

Whether  the  dominant  energy  sources  of  ULIRGs’  cores  are 
highly  obscured  normal  starbursts  with  mixed  dust/source  geom¬ 
etry  (Fig.  1  a)  or  buried  AGNs  with  centrally  concentrated  energy 
source  geometry  (Figs,  lb  and  lc)  can  be  distinguished  using  the 
equivalent  width  of  the  PAH  emission.  Since  the  PAH  equivalent 
width  (EWPah)  must  always  be  large  in  a  normal  starburst  regard¬ 
less  of  the  amount  of  dust  extinction,  a  small  EWPah  value  re¬ 
quires  a  contribution  from  a  PAH-free  continuum-emitting  energy 
source  (§  1). 

Figures  9a,  9b,  and  9c  compare  the  equivalent  widths  of  the 

6.2,  7.7,  and  11.3  pm  PAH  emission  features.  Positive  correla¬ 
tions  are  found  in  all  plots,  suggesting  that  ( 1 )  variation  of  the  rel¬ 
ative  PAH  emission  strength  is  modest,  with  some  level  of  scatter 
and  (2)  PAH-emitting  starbursts  usually  show  strong  PAH  features 
in  all  the  emission  modes  at  6.2,  7.7,  and  11.3  pm. 
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Fig.  5. — Spectra  of  all  ULIRGs  at  20bs  =  5.2—14.5  /xm,  to  zoom  in  on  the  PAH  emission  features  and  display  them  in  detail.  The  abscissa  and  ordinate  are, 
respectively,  observed  wavelength  in  /xm  and  flux  in  Jy,  both  shown  at  linear  scale.  For  IRAS  10190+1322,  the  spectra  of  the  eastern  and  western  nuclei  are  shown 
separately.  The  expected  wavelengths  of  the  6.2, 7.7,  and  11.3  /xm  PAH  emission  features,  as  well  as  6.0  /xm  H20  ice  absorption  feature,  are  indicated  as  down  arrows.  A 
question  mark  is  added  when  detection  is  unclear.  The  solid  lines  are  adopted  continuum  levels  to  estimate  the  strengths  of  the  detected  PAH  emission  features. 


Brandi  et  al.  (2006)  have  investigated  the  6.2  and  11.3  /im 
PAH  emission  equivalent  widths,  based  on  similar  definitions 
of  PAH  strength  to  ours,  in  nearby,  normal  starburst  galaxies 
without  obvious  Seyfert  signatures  in  the  optical  (18  sources), 
and  found  that  EW6.2pah  =  561  ±  73  nm  (median  553  nm)  and 
EWh  jpah  =  617  ±  187  nm  (median  597  nm).  The  scatter  is 
larger  for  the  11.3  /im  PAH  emission  feature  than  the  6.2  /jm  PAH 
feature.  No  starburst  galaxies  show  values  smaller  than  EW6.2pah  = 
459  nm  and  EWh  jpah  =  316  nm.  Since  the  PAH  equivalent 


widths  are  independent  of  the  infrared  luminosities  of  the  star- 
burst  galaxies  in  the  range  Xir  =  1 09  7— 1 01 1  6  X0  (Brandi  et  al. 
2006),  we  apply  these  values  to  our  UL1RG  sample  (Xjr  >  1O12Z0). 

The  median  values  for  the  6.2  fim  PAH  emission  equivalent 
widths  (rest-frame)  are  EW6.2pah  =  230  nm  for  the  whole  of 
our  UL1RG  sample  (LINER  +  H  n  region),  225  nm  for  LINER 
ULIRGs,  and  300  nm  for  H  n  region  ULIRGs.  The  median  11.3  /im 
PAH  equivalent  widths  are  EW i  i.3pah  =  390  nm  (LINER  +  H  n 
region),  350  nm  (LINER),  and  460  nm  (H  n  region).  The  median 
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Fig.  5 — Continued 


equivalent  widths  of  both  PAH  features  are  smaller  in  ULIRGs 
than  in  starburst  galaxies  (Brandi  et  al.  2006),  suggesting  that 
PAH-free  continuum  emission,  originating  in  AGN  -heated  hot 
dust,  contributes  more  strongly  to  the  observed  fluxes  in  ULIRGs. 
We  assume  that  if  a  PAH  equivalent  width  is  more  than  a  factor 
of~3  smaller  than  the  median  value  in  normal  starburst  galaxies, 
then  an  AGN -originating  PAH-free  continuum  contributes  im¬ 
portantly  to  the  observed  5-12  /mi  flux  and  dilutes  the  PAH  emis¬ 
sion.  With  this  criterion,  ULIRGs  with  EW6.2pah  <180  nm  and 
EW i  i.3pah  <  200  nm  are  classified  as  sources  displaying  clear 
signatures  of  luminous  AGNs.  Since  these  adopted  values  are 
clearly  below  the  minimum  values  observed  in  normal  starburst 


galaxies,  we  believe  the  criterion  is  reasonable  and  conservative 
to  pick  up  luminous  buried  AGN  candidates. 

For  the  7.7  /im  PAH  equivalent  widths,  based  on  our  defini¬ 
tion  of  the  PAH  strength  (§  4.1),  the  median  7.7  /im  PAH  equiv¬ 
alent  widths  in  our  UL1RG  sample  are  EW7.7PAH  =  490  nm 
(LINER)  and  695  nm  (H  11  region).  In  our  UL1RG  sample,  the 
median  PAH  equivalent  widths  are  smaller  in  LINER  ULIRGs 
than  in  H  11  region  ULIRGs  for  all  of  the  6.2,  7.7,  and  11.3  /im 
features.  We  classify  ULIRGs  with  EW7 7PAH  <  230  nm,  one- 
third  of  the  median  value  for  H 11  region  ULIRGs,  as  those  show¬ 
ing  luminous  buried  AGN  signatures.  Since  some  fraction  of 
H  11  region  ULIRGs  may  contain  luminous  buried  AGNs,  this 
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Fig.  5  —  Continued 


E  W7.7PAH  threshold  is  even  more  conservative  than  those  for  the 
6.2  and  11.3  //m  PAH  emission  features,  as  determined  above. 

Table  6  (cols.  [2]— [4])  presents  detection  or  nondetection 
of  buried  AGN  signatures  based  on  the  PAH  equivalent  width 
threshold.  The  fraction  of  ULIRGs  with  clearly  detectable  buried 
AGN  signatures  is  16/48  (33%;  LINER  +  H  11  region),  10/28 
(36%;  LINER)  and  6/20  (30%;  H  11  region)  based  on  observed 
EW6.2PAH,  or  20/48  (42%;  LINER  +  H  11  region),  13/28  (46%; 
LINER),  and  7/20  (35%;  H  11  region)  based  on  EW6.2pah  after 
correction  for  the  ice  and  HAC  absorption  feature.  The  fraction 
is  4/48  (8%;  LINER  +  H 11  region),  3/28  (11%;  LINER),  and  1/20 
(5%;  H  11  region)  based  on  observed  EW7.7PAH,  12/48  (25%; 


LINER  +  H  11  region),  8/28  (29%;  LINER),  and  4/20  (20%;  H  11 
region)  based  on  observed  EWhjpah-  The  fraction  of  ULIRGs 
where  very  low  PAH  equivalent  widths  are  found  in  any  of  the 
features  is  19/48  (40%;  LINER  +  Hu  region),  13/28  (46%;  LINER), 
and  6/20  (30%;  H 11  region).  If  the  corrected  EW6  2pah  values  are 
adopted,  it  is  22/48  (46%;  LINER  +  H  11  region),  15/28  (54%; 
LINER),  and  7/20  (35%;  H  11  region). 

Since  we  adopted  a  conservative  threshold,  this  method  pro¬ 
vides  only  candidate  ULIRGs  where  PAH-free  continuum  emis¬ 
sion  from  buried  AGNs  contribute  very  importantly  to  the  observed 
IRS  spectra.  In  starburst/buried-AGN  composite  ULIRGs  (Fig.  lc), 
PAH  emission  from  starbursts  surrounding  the  central  buried 
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AGNs  is  less  obscured  and  less  flux-attenuated,  making  the  ob¬ 
served  PAH  equivalent  widths  relatively  large,  particularly  when 
AGNs  are  deeply  buried.  Such  buried  AGNs  should  be  searched 
for  using  other  methods. 

5.2.2.  Absolute  Optical  Depths  of  Dust  Absorption  Features 

In  a  normal  starburst  with  mixed  dust/source  geometry  (Fig.  la), 
the  foreground,  less-obscured,  less-attenuated  emission  (which 
shows  only  weak  dust  absorption  features)  dominates  the  ob¬ 
served  flux,  and  therefore  the  observed  optical  depths  of  dust 
absorption  features  cannot  exceed  a  certain  threshold,  unless 
very  unusual  dust  composition  patterns  are  assumed  (Imanishi 
&  Maloney  2003;  Imanishi  et  al.  2006a).  In  a  buried  AGN,  with 
centrally  concentrated  source  geometry  (Fig.  lb),  the  so-called 
foreground  screen  dust  model  is  applicable,  and  the  observed 
optical  depths  can  be  arbitrarily  large.  Hence,  detection  of  strong 
dust  absorption  features,  with  optical  depths  substantially  larger 
than  the  upper  limit  set  by  the  mixed  dust/source  geometry,  ar¬ 
gues  for  a  foreground  screen  dust  geometry,  as  expected  from  a 
buried  AGN  (Imanishi  &  Maloney  2003;  Imanishi  et  al.  2006a). 

To  differentiate  the  geometry  using  the  optical  depths  of  dust 
absorption  features,  the  9.7  pm  silicate  dust  absorption  is  par¬ 
ticularly  useful,  because  (1)  it  is  much  stronger  than  the  18  pm 
silicate  feature  (Fig.  6a)  and  (2)  its  intrinsic  strength  has  been 
widely  investigated  in  Galactic  obscured  objects  ( Roche  &  Aitken 
1984, 1985).  Assuming  the  Galactic  dust  extinction  curve,  the  up¬ 
per  limit  for  the  optical  depth  of  the  9.7  pm  silicate  absorption 
feature  in  the  mixed  dust/source  geometry  can  be  derived.  The 
Appendix  presents  detailed  calculations  and  we  obtain  a  maxi¬ 
mum  value  of  T9  .7  <  1  .7  for  mixed  dust/source  geometry.  In  fact, 
none  of  the  starburst  galaxies  (ZIR  <  101L6  Z0)  lacking  optical 
Seyfert  signatures  so  far  observed  with  Spitzer  IRS  (Brandi  et  al. 
2006)  show  values  that  exceed  the  threshold  (r^  7  =  1.7).  Thus, 
it  is  reasonable  to  assume  that  ULIRGs  with  r'g  7  substantially 
larger  than  1 .7  possess  buried  AGNs  with  centrally  concentrated 
energy  source  geometry. 

In  ULIRGs,  it  is  possible  that  the  dust  extinction  curve  is  very 
different  from  the  Galactic  one,  because  dust  coagulation  in  the 
high-density  environment  in  ULIRGs’  cores  could  increase  typi¬ 
cal  dust  grain  size  up  to  a  few  microns  (Laor  &  Draine  1993; 
Maiolino  et  al.  2001a,  2001b;  Imanishi  2001).  However,  we  be¬ 


lieve  this  possibility  has  minimal  effect  on  our  conclusions  about 
the  geometry,  as  long  as  we  adopt  Tg  7  (peak  optical  depth  of  the 
9.7  pm  feature,  relative  to  the  nearby  continuum  just  outside 
the  feature;  see  §  4.2),  rather  than  Tgj  (peak  optical  depth  of  the 
9.7  pm  feature  and  continuum,  relative  to  continuum  outside  8- 
30  /im;  see  §  4.2).  Figure  6b  shows  the  absorption  optical  depth 
of  astronomical  silicate  dust  grains  of  three  different  sizes.  The 
absolute  absorption  optical  depth  at  the  peak  of  the  9.7  pm  fea¬ 
ture  becomes  stronger  with  increasing  dust  grain  size.  However, 
both  the  absorption  optical  depth  inside  the  9.7  pm  feature  and 
the  nearby  continuum  just  outside  the  feature  increase  in  a  very 
similar  way  with  increasing  dust  size,  so  that  the  strength  of  9.7  pm 
feature  relative  to  the  nearby  continuum  (Tg  7),  is  virtually  un¬ 
changed  with  varying  dust  grain  size.  The  argument  of  the  geom¬ 
etry  based  on  r9'  7  values  are  therefore  robust  to  possible  increase 
of  dust  grain  size  in  ULIRGs’  cores.  We  conservatively  take 
ULIRGs  with  Tg  7  >  2  as  candidates  for  harboring  luminous  cen¬ 
trally  concentrated  buried  AGNs. 

Individual  ULIRGs  with  clearly  detectable  buried  AGN  sig¬ 
natures  based  on  large  r9'  7  values  (>2;  Table  4)  are  marked  with 
open  circles  in  Table  6  (col.  [5]).  The  median  T9  7  values  are  1.8 
(LINER +  H 11  region),  2.1  (LINER),  and  1.6  (Hii  region).  The 
fraction  of  ULIRGs  with  t9.7  >  2,  indicative  of  buried  AGNs  with 
centrally  concentrated  energy  source  geometry,  is  20/48  (42%; 
LINER  +  H  11  region),  14/28  (50%;  LINER),  and  6/20  (30%;  H  11 
region).  The  fraction  of  ULIRGs  showing  buried  AGN  signatures 
is  again  higher  in  LINER  ULIRGs  than  H  11  region  ULIRGs,  as 
was  previously  found  based  on  the  PAH-equivalent- width  method 
(§  5.2.1). 

We  must  raise  three  cautions  regarding  this  method.  First,  a 
T9  7  value  smaller  than  1.7  does  not  necessarily  preclude  the  pres¬ 
ence  of  a  luminous  AGN,  because  a  weakly  obscured  AGN,  of 
course,  shows  a  small  T9  7  value.  The  best  example  is  IRAS 
01004—2237  (t9.7  ~  0.4),  which  shows  very  low  PAH  emis¬ 
sion  equivalent  widths,  as  expected  from  an  AGN-dominated 
source. 

Next,  we  used  the  threshold  of  Tg  .7  >  2,  based  on  the  observed 
values.  For  ULIRGs  with  detectable  PAH  emission  from  mod¬ 
estly  obscured  starbursts,  the  observed  fluxes  are  a  superposition 
of  the  emission  from  the  highly  obscured  cores  and  the  starburst 
emission  (Fig.  lc),  and  so  the  observed  Tg7  values  are  reduced 


TABLE  3 
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U) 

EW6.2PAH 

(nm) 

(2) 

EW7  7PA[,a 
(nm) 

(3) 

EWn  3PAH 
(nm) 

(4) 

EW33pah 

(nm) 

(5) 

E6.2PAH 

(1042  ergs  s_1) 
(6) 

T  a 

^7.7PAH 

(1042  ergs  s-1) 
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Eu.3pah 
(1042  ergs  s-1) 
(8) 

E6.2PAH/£lR 

(10-3) 

(9) 

Ell.3PAn/ElR 

(10-3) 

(10) 

IRAS  00188-0856 . 

85  (~40) 

305 

290 

50 

3.5 

17 

3.4 

0.4 

0.4 

IRAS  00482-2721 . 

260 

710 

465 

1.8 

4.8 

0.9 

0.5 

0.2 

IRAS  03250+1606 . 

325 

740 

570 

80 

7.0 

18 

4.6 

1.6 

1.1 

IRAS  04103-2838 . 

165 

365 

120 

8.4 

21 

6.6 

1.4 

1.1 

IRAS  08572+3915 . 

<15 

<100 

<30 

<5 

<1.9 

<18 

<0.5 

<0.4 

<0.1 

IRAS  09039+0503 . 

.  325  (—200) 

625 

580 

95 

3.7 

11 

2.2 

0.8 

0.5 

IRAS  09116+0334 . 

360 

730 

630 

75 

9.4 

27 

7.4 

1.9 

1.5 

IRAS  09539+0857 . 

145 

405 

390 

65 

2.5 

12 

1.3 

0.6 

0.3 

IRAS  10378+1108 . 

110  (—50) 

450 

210 

40 

2.2 

14 

1.5 

0.3 

0.2 

IRAS  10485-1447 . 

.  225  (—150) 

400 

150 

80 

3.7 

11 

1.1 

0.6 

0.2 

IRAS  10494+4424 . 

.  345  (—250) 

715 

495 

110 

6.1 

19 

3.1 

1.1 

0.6 

IRAS  11095-0238 . 

20  (—10) 

250 

110 

150 

1.0 

17 

0.9 

0.2 

0.1 

IRAS  11130-2659 . 

.  95  (—60) 

490 

280 

1.9 

13 

1.0 

0.4 

0.2 

IRAS  12112+0305 . 

390 

565 

425 

100+160 

6.1 

15 

3.2 

0.8 

0.4 

IRAS  12127-1412 . 

5  (~3) 

95 

20 

0 

0.7 

11 

0.6 

0.1 

0.1 

IRAS  12359-0725 . 

220 

445 

165 

75 

3.1 

8.8 

1.1 

0.6 

0.2 

IRAS  13335-2612 . 

310 

700 

610 

6.6 

17 

5.1 

1.5 

1.2 

IRAS  14252-1550 . 

225 

685 

320 

60+60 

3.6 

12 

2.3 

0.6 

0.4 

IRAS  14348-1447 . 

.  260  (—110) 

540 

500 

70+110 

5.9 

19 

3.6 

0.8 

0.5 

Arp  220 . 

.  215  (—120) 

410 

315 

80 

1.5 

5.1 

0.8 

0.3 

0.2 

IRAS  16090-0139 . 

110  (—70) 

265 

325 

75 

8.1 

30 

5.6 

0.7 

0.5 

IRAS  16468+5200 . 

230 

470 

435 

120 

2.9 

11 

2.0 

0.7 

0.5 

IRAS  16487+5447 . 

250 

580 

570 

70 

2.8 

9.6 

2.2 

0.5 

0.4 

IRAS  17028+5817 . 

375 

800 

595 

120 

4.6 

14 

2.8 

0.9 

0.6 

IRAS  17044+6720 . 

55 

165 

30 

10 

3.9 

13 

1.0 

0.8 

0.2 

IRAS  21329-2346 . 

230 

635 

370 

50 

3.0 

12 

1.3 

0.6 

0.3 

IRAS  23234+0946 . 

220 

695 

445 

75 

3.1 

11 

2.0 

0.7 

0.4 

IRAS  23327+2913 . 

255 

485 

175 

45 

2.5 

6.9 

1.4 

0.6 

0.3 

IRAS  00091-0738 . 

50  (—30) 

350 

255 

1.3 

15 

1.5 

0.2 

0.3 

IRAS  00456-2904 . 

385 

750 

545 

11 

27 

7.6 

1.9 

1.4 

IRAS  01004-2237 . 

45 

190 

8 

3.1 

19 

0.5 

0.4 

0.1 

IRAS  01166-0844 . 

30 

430 

190 

0.5 

8.5 

0.7 

0.1 

0.2 

IRAS  01298-0744 . 

10  (—6) 

380 

155 

0.7 

37 

1.1 

0.1 

0.1 

IRAS  01569-2939 . 

115 

555 

255 

3.1 

20 

2.0 

0.5 

0.3 

IRAS  02411+0353 . 

300 

695 

375 

13 

38 

7.9 

2.1 

1.3 

IRAS  10190+1322E . 

375 

680 

510 

95 

4.4 

12 

3.0 

1.2 

0.8 

IRAS  10190+1322W . 

375 

745 

685 

50 

4.6 

11 

3.6 

1.2 

0.9 

IRAS  11387+4116 . 

265 

680 

475 

75 

4.2 

15 

3.9 

0.7 

0.7 

IRAS  11506+1331  . 

195  (—120) 

540 

305 

95 

8.5 

26 

3.4 

1.1 

0.4 

IRAS  13509+0442 . 

490 

795 

590 

135 

9.7 

22 

4.7 

1.5 

0.7 

IRAS  13539+2920 . 

355 

755 

600 

85 

7.6 

22 

4.4 

1.9 

1.1 

IRAS  14060+2919 . 

340 

860 

530 

150 

8.6 

27 

6.5 

2.0 

1.5 

IRAS  15206+3342 . 

225 

525 

230 

55 

12 

32 

7.7 

2.0 

1.3 

IRAS  15225+2350 . 

145  (—100) 

290 

130 

40 

5.0 

16 

2.0 

1.0 

0.4 

IRAS  16474+3430 . 

390 

820 

575 

105 

9.2 

27 

5.8 

1.7 

1.1 

IRAS  20414-1651 . 

330 

705 

465 

75 

3.2 

10 

1.7 

0.5 

0.3 

IRAS  21208-0519 . 

505 

810 

420 

100 

6.0 

14 

2.7 

1.5 

0.7 

IRAS  22206-2715 . 

210 

975 

565 

3.1 

13 

2.7 

0.6 

0.5 

IRAS  22491-1808 . 

325 

695 

460 

3.9 

12 

2.6 

0.8 

0.5 

Notes. — Col.  (1):  Object  name.  Col.  (2):  Rest-frame  equivalent  width  of  the  6.2  fim  PAH  emission.  The  values  in  parentheses  are  approximate  ones  after  correction 
for  the  ice  and  HAC  absorption  features.  Col.  (3):  Rest-frame  equivalent  width  of  the  7.7  fim  PAH  emission.  Col.  (4):  Rest-frame  equivalent  width  of  the  1 1 .3  fim  PAH 
emission.  Col.  (5):  Rest-frame  equivalent  width  of  the  3.3  /m l  PAH  emission,  taken  from  Imanishi  et  al.  (2006a).  Ellipses  signify  no  data.  Col.  (6):  Luminosity  of  the 
6.2  fim  PAH  emission  in  units  of  1042  ergs  s-1.  Col.  (7):  Luminosity  of  the  7.7  fim  PAH  emission  in  units  of  1042  ergs  s-1.  Col.  (8):  Luminosity  of  the  1 1 .3  fim  PAH  emission 
in  units  of  1042  ergs  s”1.  Col.  (9):  The  6.2  fim  PAH-to-infrared  luminosity  ratio  in  units  of  10-3.  The  ratio  for  normal  starbursts  with  modest  dust  obscuration 
(Ay  <  20—30  mag)  is  ~3.4  x  10-3  (Peeters  et  al.  2004).  Col.  (10):  The  1 1 .3  fim  PAH-to-infrared  luminosity  ratio  in  units  of  10-3.  The  ratio  for  normal  starbursts  with 
modest  dust  obscuration  (Ay  <  20—30  mag)  is  ~1.4x  10-3  (Soifer  et  al.  2002). 

a  We  regard  flux  excess  at  2rest  =  7.3— 8.1  fim  above  an  adopted  continuum  level  as  7.7  fim  PAH  emission,  to  reduce  the  effects  of  the  strong  9.7  fim  silicate  dust 
absorption  feature.  Our  definition  is  different  from  those  in  many  previous  literatures. 
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Fig.  6. — Silicate  dust  absorption  profile,  (a)  Solid  lines  are  the  observed  silicate  absorption  profiles  of  the  Galactic  stars  GCS  3-1  and  WR  98a  (Chiar  &  Tielens  2006). 
The  dashed  line  is  the  theoretical  absorption  profile  of  “astronomical  silicate”  grains  with  0.1  //m  in  size  (Draine  &  Lee  1984).  All  profiles  are  normalized  at  the  9.7  /.tm 
absorption  peak.  The  silicate  absorption  profile  extends  from  8  to  >30  pm.  (b)  Theoretical  absorption  profile  of  “astronomical  silicate”  grains  (Draine  &  Lee  1984)  of  dif¬ 
ferent  sizes,  1  //m  ( dotted  line),  0.1  //in  (dashed  line),  and  0.01  //m  ( solid  line).  Adopted  from  http://www.astro.princeton.edu/~draine/dust/dust.diel.html.  The  ordinate  is  in 
arbitrary  units  but  is  proportional  to  the  actual  optical  depth.  In  this  figure,  only  absorption  is  taken  into  account,  and  scattering  is  neglected,  because  scattering  is  negligible  in 
this  wavelength  range  (see  the  above  Web  page). 


compared  to  the  Tg  7  values  only  toward  the  highly  obscured 
energy  sources  at  ULIRGs  ’  cores.  Since  subtraction  of  the  mod¬ 
estly  obscured  starburst  component,  using  a  starburst  template 
spectrum,  introduces  an  additional  uncertainty,  we  do  not  at¬ 
tempt  it.  Thus,  we  may  miss  ULIRGs  that  show  large  r9'  7  values 
(>2)  toward  the  cores  but  small  observed  r9'  7  values  (<2)  due  to 
the  starburst  dilution. 

Finally,  a  starburst  nucleus  obscured  by  foreground  absorbing 
dust  (Fig.  1  d)  can  also  show  a  r9  7  value  larger  than  the  threshold 
determined  by  the  mixed  dust/source  geometry,  because  in  this 
case,  a  foreground  screen  dust  model  is  applicable.  Assuming 
Galactic  dust  properties,  r97  >  2  corresponds  to  Ay  >  20-40  mag 
(Roche  &  Aitken  1984,  1985).  If  ULIRGs  are  the  mergers  of  gas- 
rich  spiral  galaxies  (Sanders  et  al.  1 988a),  gas  and  dust  in  the  orig¬ 
inal  galaxies  can  be  quickly  transferred  to  the  nuclear  regions 
(Barnes  &  Hemquist  1996;  Mihos  &  Hemquist  1996).  The  bulk 
of  possibly  newborn  stars  and  resulting  dust  are  also  expected  to 
distribute  in  the  nuclear  regions.  In  fact,  most  of  high-density  gas, 
dust,  and  energy  sources  (starbursts  and/or  AGNs)  in  ULIRGs 
are  observed  to  be  concentrated  to  nuclear  cores  with  <kpc  scale 
(Downes  &  Solomon  1998;  Soifer  et  al.  2000;  Imanishi  et  al. 
2006b).  The  amount  of  gas/dust  in  ULIRGs’  hosts  at  greater  than 
kpc  scale  (which  are  morphologically  disturbed)  is  unlikely  to  sub¬ 
stantially  exceed  that  of  the  original  merging  gas/dust-rich  galax¬ 
ies.  In  our  Galaxy,  a  gas-rich  spiral,  the  A  v  value  is  estimated  to  be 
>30  mag  in  an  equatorial  direction  (Rieke  et  al.  1989)  and  should 
be  much  smaller  in  other  directions.  It  is  thus  not  easy  to  produce 
such  a  large  A  y  value  (>20-40  mag)  with  dust  in  a  host  galaxy, 
unless  the  host  is  viewed  from  a  direction  with  particularly  large 
dust  column  density.  Although  a  small  fraction  of  ULIRGs  with 
r9  7  >2  may  have  this  geometry  (Fig.  Id),  it  is  very  unlikely  that 
the  majority  do. 

5.2.3.  Strong  Dust  Temperature  Gradients 

A  buried  AGN  with  centrally  concentrated  energy  source 
geometry  (Fig.  lb)  should  show  a  strong  dust  temperature  gra¬ 
dient  (Fig.  2),  whereas  a  normal  starburst  nucleus  with  mixed 
dust/source  geometry  does  not,  whether  it  is  obscured  by  fore¬ 
ground  dust  (Fig.  \d)  or  not  (Fig.  In).  The  presence  of  a  strong 
dust  temperature  gradient  thus  supports  the  buried  AGN  scenario. 


and  it  can  be  detected  by  comparing  the  optical  depths  of  dust 
absorption  features  at  different  infrared  wavelengths  (Dudley 
&  Wynn- Williams  1997;  Imanishi  2000,  2006a;  Imanishi  & 
Maloney  2003).  As  Figure  2  illustrates,  in  a  buried  AGN,  the 
dust  column  density  estimated  from  observations  at  2rest  >  3  pm 
probes  deeper  inside  of  the  obscuring  material  as  the  wavelength 
is  decreased  [i.e.,  A  y(3  pm)  >  Ay(  10  pm)  >  Av(  20  /.tm)].  In  this 
geometry,  dust  with  a  fixed  temperature  range  from  a  fixed  layer 
contributes  strongly  to  the  observed  continuum  emission  at  each 
wavelength,  with  a  small  contribution  from  other  components 
(Fig.  2).  Thus,  in  a  starburst/buried  AGN  composite  (Fig.  lc), 
PAH-free  continuum  emission  from  the  AGN-heated  dust  should 
be  detectable  even  if  its  optical  depth  is  larger  than  unity  at  a  par¬ 
ticular  observed  wavelength,  provided  the  buried  AGN  is  intrin¬ 
sically  sufficiently  luminous  that  the  flux-attenuated  (e-r;  r  >  1) 
AGN-heated  dust  emission  contributes  significantly  to  the  ob¬ 
served  flux.  This  has  been  demonstrated  in  the  3-4  pm  spectra  of 
the  ULIRGs  IRAS  08572+3915  and  UGC  5101  (t3_4  ;,m  >  5; 
Imanishi  2000,  2001;  Imanishi  &  Maloney  2003). 

On  the  other  hand,  in  a  normal  starburst  (Fig.  1  a),  no  sig¬ 
nificant  dust  temperature  gradient  is  expected,  and  so  the  esti¬ 
mated  dust  extinction  should  be  similar  at  different  wavelengths 
[Av( 3  pm)  ~  Ay(10  pm)  ~  Av(20  /tm)].  Since  in  a  normal  star- 
burst,  emission  arises  uniformly  within  an  emitting  volume,  we 
can  typically  probe  emission  with  optical  depth  less  than  about 
unity.  If  the  extinction  is  large  enough  that  the  short  wavelengths 
do  not  probe  the  whole  emitting  region,  the  estimated  dust  extinc¬ 
tion  would  decrease  with  decreasing  wavelength  [Ay (3  pm)  < 
Av(  10  pm)  <  Av(  20  /im)].  Hence,  the  relations  among  A  K(  3  pm), 
Ay(  10  pm),  and  Ay (20  pm)  are  totally  different  between  a  bur¬ 
ied  AGN  (Fig.  lb)  and  a  normal  starburst  (Fig.  la). 

For  a  buried  AGN,  the  t{s/tq  7  or  Tj'g/r9  7  ratio  in  the  Spitzer 
IRS  spectrum  becomes  small  compared  to  the  value  in  the  absence 
of  a  strong  dust  temperature  gradient.  The  Tj'8Ar9  7  (tiVt9.7)  rat*° 
is  also  insensitive  to  possible  changes  in  dust  grain  size  in 
ULIRGs,  because  the  relative  strength  of  r(8  (t"s)  and  r9  7  varies 
little  with  varying  dust  grain  size  (Fig.  6b).  Furthermore,  since 
both  absorption  features  are  due  to  silicate  dust  grains,  the  ratio 
is  unaffected  by  possibly  different  dust  composition  patterns  in 
individual  galaxies. 


TABLE  4 

Optical  Depth  of  the  9.7  fim  and  18  fim  Silicate  Dust  Absorption  Features 


Object 

(1) 

T9.7 

(2) 

r9.7 

(3) 

t!s 

(4) 

<g 

(5) 

T<  lT' 

'18 ''9.7 

(6) 

rlVr9.7 

(7) 

IRAS  00188-0856 . 

.  3.1 

2.5 

0.55 

0.65 

0.22 

0.26 

IRAS  00482-2721  . 

.  2.6 

2.1 

0.6 

0.75 

0.29 

0.36 

IRAS  03250+1606 . 

.  1.8 

1.6 

IRAS  04103-2838 . 

.  0.5 

0.6 

IRAS  08572+3915 . 

.  4.1 

3.8 

0.9 

1.0 

0.24 

0.26 

IRAS  09039+0503 . 

.  2.6 

2.2 

0.9 

1.0 

0.41 

0.45 

IRAS  09116+0334 . 

.  1.5 

1.1 

IRAS  09539+0857 . 

.  3.8 

3.5 

1.3 

1.6 

0.37 

0.46 

IRAS  10378+1108 . 

.  2.6 

2.4 

0.45 

0.55 

0.19 

0.23 

IRAS  10485-1447 . 

.  3.0 

3.0 

0.8 

0.95 

0.27 

0.32 

IRAS  10494+4424 . 

.  2.2 

1.7 

IRAS  11095-0238 . 

.  3.5 

3.4 

0.8 

0.9 

0.23 

0.26 

IRAS  11130-2659 . 

.  3.5 

3.1 

1.0 

1.2 

0.32 

0.39 

IRAS  12112+0305 . 

.  2.3 

1.7 

IRAS  12127-1412 . 

.  3.0 

2.5 

0.6 

0.7 

0.24 

0.28 

IRAS  12359-0725 . 

.  2.2 

1.8 

0.5 

0.55 

0.28 

0.31 

IRAS  13335-2612 . 

.  1.6 

1.4 

IRAS  14252-1550 . 

.  1.6 

1.4 

IRAS  14348-1447 . 

.  2.8 

2.0 

Ar)i  220 . 

.  3.7 

3.2 

0.95 

1.1 

0.30 

0.34 

IRAS  16090-0139 . 

.  3.0 

2.6 

0.7 

0.8 

0.27 

0.31 

IRAS  16468+5200 . 

.  3.0 

2.5 

0.65 

0.8 

0.26 

0.32 

IRAS  16487+5447 . 

.  2.4 

1.8 

IRAS  17028+5817 . 

.  2.0 

1.5 

IRAS  17044+6720 . 

.  1.7 

1.8 

0.4 

0.5 

0.22 

0.28 

IRAS  21329-2346 . 

.  2.7 

2.2 

IRAS  23234+0946 . 

.  2.2 

1.8 

IRAS  23327+2913 . 

.  1.5 

1.3 

IRAS  00091-0738 . 

.  3.2 

3.2 

1.05 

1.25 

0.33 

0.39 

IRAS  00456-2904 . 

.  1.5 

1.2 

IRAS  01004-2237 . 

.  0.4 

0.4 

IRAS  01166-0844 . 

.  3.2 

3.0 

0.75 

0.9 

0.25 

0.30 

IRAS  01298-0744 . 

.  4.3 

4.0 

1.1 

1.35 

0.27 

0.34 

IRAS  01569-2939 . 

.  3.2 

2.8 

1.0 

1.2 

0.36 

0.43 

IRAS  02411+0353 . 

.  1.3 

1.0 

IRAS  10190+1322“ . 

.  1.7 

1.3 

IRAS  11387+4116 . 

.  1.6 

1.1 

IRAS  11506+1331 . 

.  2.6 

2.3 

0.7 

0.85 

0.30 

0.37 

IRAS  13509+0442 . 

.  2.0 

1.6 

IRAS  13539+2920 . 

.  1.8 

1.6 

IRAS  14060+2919 . 

.  1.5 

1.0 

IRAS  15206+3342 . 

.  0.9 

0.6 

IRAS  15225+2350 . 

.  2.5 

2.2 

0.65 

0.75 

0.29 

0.34 

IRAS  16474+3430 . 

.  1.8 

1.5 

IRAS  20414-1651 . 

.  2.4 

1.9 

IRAS  21208-0519 . 

.  1.5 

1.4 

IRAS  22206-2715 . 

.  2.4 

1.7 

IRAS  22491-1808 . 

.  2.0 

1.5 

Notes. — Col.  (1):  Object  name.  Col.  (2):  r9  7  is  an  optical  depth  of  the  9.7  fim  silicate  dust  absorption  feature,  against  a 
power-law  continuum  determined  from  data  points  at  2rest  =5.6  and  7. 1  fim  and  20bS  =  34—35  /im,  shown  as  dashed  lines 
in  Figs.  3  and  4.  Col.  (3):  7  is  an  optical  depth  of  the  9.7  /im  silicate  dust  absorption  feature,  against  a  power-law  con¬ 

tinuum  determined  from  data  points  at  2rest  =  7.1  and  14.2  /im,  shown  as  dotted  lines  in  Figs.  3  and  4.  Once  the  continuum 
levels  are  fixed,  the  uncertainty  of  r9  7  is  <5%  for  ULIRGs  with  large  r9  7  values  (>2)  and  can  be  ~10%  for  ULIRGs  with 
small  r97.  Col.  (4):  Tj8  is  an  optical  depth  of  the  18  /im  silicate  dust  absorption  feature,  against  a  power-law  continuum 
determined  from  data  points  at  2^  —  14.2  and  24  /im,  shown  as  dotted  lines  in  Figs.  3  and  4.  Once  the  continuum  is  fixed,  the 
statistic  uncertainty  of  r18  is  <10%,  because  the  value  is  estimated  only  for  ULIRGs  with  clearly  detectable  18  fim  silicate 
absorption  features.  Col.  (5):  t[ g  is  an  optical  depth  of  the  1 8  /im  silicate  dust  absorption  feature,  against  a  power-law  continuum 
determined  from  data  points  at  2rest  =  14.2  and  29  /im,  shown  as  dash-dotted  lines  in  Figs.  3  and  4.  Once  the  continuum  is  fixed, 
the  statistic  uncertainty  is  <  10%.  Col.  (6):  r[g  /r9  7  ratio,  only  for  ULIRGs  with  clearly  detectable  18  fim  silicate  absorption.  The 
uncertainty  is  <  10%  once  the  continuum  is  fixed.  Col.  (7):  r(g/r9  7  ratio,  only  for  ULIRGs  with  clearly  detectable  18  fim  silicate 
absorption. 

a  To  estimate  the  optical  depths  of  silicate  dust  absorption  features,  SL  spectra  of  the  eastern  and  western  nuclei  are 
summed. 
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Fig.  7. — IRS  spectra  at  20bs  —  5.2—9  fim  for  ULIRGs  displaying  clear  6.0  fim  H2O  absorption  feature  (marked  with  “6.0  fim  ice  absorption”  in  the  first  15  plots). 
Spectra  of  the  next  five  ULIRGs  [marked  with  “6.0  fim  ice  absorption  (?)”]  show  possible  signatures  of  the  6.0  fim  ice  absorption  feature.  Spectra  of  the  last  two 
ULIRGs,  IRAS  17044+6720  and  00456—2904,  marked  with  “no  6.0  fim  ice  absorption”  are  shown  as  examples  with  undetected  ice  absorption  features.  The  abscissa  is 
observed  wavelength  in  fim  in  linear  scale,  and  the  ordinate  is  flux  in  Jy  in  decimal  logarithmic  scale. 


For  ULIRGs  with  large  contributions  from  modestly  obscured 
PAH-emitting  starbursts  to  observed  spectra,  the  18  fim  sili¬ 
cate  dust  absorption  features  are  diluted  by  the  PAH  complex 
at  Arest  =  16.4  and  17.12  fim,  as  well  as  H2  emission  at  Arest  = 
17.04  and  [Sin]  emission  at  Arest  =  18.71  fim  (Smith  et  al.  2004). 
Subtraction  of  a  template  spectrum  of  the  PAH-emitting  starburst 
introduces  an  additional  ambiguity  in  the  estimate  of  the  intrinsic 
Tj g  (r"g)  value  toward  the  ULIRG’s  core.  To  minimize  the  un¬ 
certainties,  we  investigate  the  T['8/t9  7  (ti'8/t9  7)  ratios  only  for 
ULIRGs  where  the  contributions  from  PAH-emitting  normal 
starbursts  are  small  and  the  18  /tm  silicate  features  are  clearly 
seen  in  absorption  in  the  observed  spectra.  The  main  aim  of  ap¬ 
plying  this  Tj'8/t9  7  (Yj'g/rg  7)  method  to  these  ULIRGs  is  to  fur¬ 
ther  strengthen  the  case  for  buried  AGNs  previously  suggested 
by  the  large  r9  7  values,  rather  than  to  discover  new  buried  AGN 
candidates. 

To  investigate  the  presence  of  a  strong  dust  temperature  gra¬ 
dient  in  a  ULIRG,  based  on  a  small  observed  (T]'g/T9  7) 

ratio,  we  first  have  to  know  its  intrinsic  ratio  in  the  absence  of  a 
temperature  gradient.  Figure  6a  presents  the  observed  profiles 


of  two  Galactic  obscured  stars,  GCS3-I  and  WR  98a  (Chiar  & 
Tielens  2006).  They  are  thought  to  be  mostly  obscured  by  fore¬ 
ground  interstellar  medium  far  away  from  the  central  illumi¬ 
nating  stars,  rather  than  circumstellar  material,  so  no  significant 
temperature  gradient  is  expected.  To  conservatively  argue  that  a 
strong  dust  temperature  gradient  is  found,  based  on  a  small  ob¬ 
served  Tj'g/rg  7  ("Tj'g/rg  7)  ratio,  we  have  adopted  the  profile  show¬ 
ing  the  smallest  intrinsic  (n'g/r^)  ratio  as  a  template. 

The  intrinsic  ratio  of  GCS3-I  is  estimated  to  be  Ti'8/t9  7  =  0.3 
and  Tjg/T97  =  0.35.  The  ratios  of  WR  98a,  and  other  Galactic 
obscured  stars  for  which  silicate  profiles  are  investigated  in 
detail  (OH-1R127. 8+0.0,  WR  1 1 8,  and  WR  1 12;  Chiar  &  Tielens 
2006),  are  larger  than  that  of  GCS3-I,  so  that  adopting  the  ob¬ 
served  GCS3-I  profile  as  a  template  is  appropriate  here.  If  an  ob¬ 
served  Tj'g  /r9  7  or  t{8/t9  7  ratio  is  substantially  smaller  than  the 
adopted  intrinsic  value  (r/g/rg  7  =  0.3  and  Tj'8/T9  7  =  0.35),  then 
the  most  reasonable  explanation  is  the  presence  of  a  strong  dust 
temperature  gradient. 

Table  4  (cols.  [6]  and  [7])  summarizes  the  'r1,g/'r9  7  and  T[8/t9  7 
ratios  for  ULIRGs  showing  clear  18  fim  silicate  absorption 
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Fig.  7 — Continued 


(mostly  Tg  7  >  2).  For  a  significant  fraction  of  such  ULIRGs,  the 
observed  ratios  are  smaller  than  the  conservatively  assumed 
intrinsic  one,  suggesting  the  presence  of  strong  dust  temperature 
gradients.  This  provides  additional  evidence  for  these  ULIRGs 
for  centrally  concentrated  buried  AGNs,  previously  suggested 
by  low  PAH  equivalent  widths  and/or  large  Tg  7  values  (>2). 

ULIRGs  showing  detectable  signatures  of  strong  dust  tem¬ 
perature  gradients  are  marked  as  open  circles  in  Table  6  (col.  [6]). 
The  fraction  of  ULIRGs  with  detectable  signatures  of  strong  dust 
temperature  gradients  is  7/48  (15%;  LINER  +  H  n  region),  6/28 
(21%;  LINER),  and  1/20  (5%;  H  n  region).  The  fraction  of  ULIRGs 
with  detectable  buried  AGN  signatures  is  again  higher  in  LINER 
ULIRGs  than  H  n  region  ULIRGs. 

Again,  there  are  several  caveats.  First,  we  adopted  the  observed 
profile  of  GCS3-I  as  the  template,  since  the  r(s/Tg  7  ratio  in  GCS3-I 
is  the  lowest  among  five  Galactic  sources  so  far  investigated  in 
detail.  A  source  with  an  even  smaller  intrinsic  +1V+9.  7  ratio  may 
be  discovered  among  Galactic  sources  in  future  observations.  If 
ULIRGs  showed  intrinsically  smaller  r'^/Tgj  ratios  than  the 
adopted  ratio  from  GCS3-I,  then  the  small  observed  ratios  in 
ULIRGs  would  not  necessarily  require  the  strong  dust  tempera¬ 
ture  gradients.  However,  for  IRAS  08572+3915,  a  bright  UL1RG 


with  a  small  Tj'8 /Tg  7  ratio  (Table  4),  Spoon  et  al.  (2006)  indepen¬ 
dently  found  the  signature  of  a  strong  dust  temperature  gradient, 
based  on  the  fact  that  the  observed  optical  depths  of  crystalline 
silicate  dust  absorption  at  Lrest  ~  11,  16,  18.5,  23,  and  28  pm  are 
systematically  weaker  at  the  longer  wavelengths  than  expected 
(their  Fig.  2,  bottom).  A  similar  trend  was  found  also  in  the  aver¬ 
age  spectrum  of  ULIRGs  with  Tg1>  2  (Spoon  et  al.  2006,  their 
Fig.  5).  We  thus  believe  that  the  small  rfg/rg  7  ratios  observed  in 
several  ULIRGs  are  most  naturally  explained  by  strong  dust  tem¬ 
perature  gradients. 

Next,  although  we  investigated  the  Ti8/t9.7  ratios  only  for 
ULIRGs  with  small  PAH  contributions  from  modestly  obscured 
starbursts,  these  small  contributions  could  reduce  both  +18  and 
Tgj  (§  5.2.2),  and  slightly  alter  the  r[s/Tg  7  ratios,  compared  to  a 
pure  spectrum  of  the  highly  obscured  energy  sources  in  the  cores. 
However,  the  average  spectrum  of  starburst  galaxies  shows  the 
30-5  pm  (rest-frame)  flux  ratio  to  be  ^900  (Brandi  et  al.  2006), 
while  the  majority  of  the  ULIRGs  with  small  r/g/rg  7  show  flux 
ratios  larger  than  >900  (Figs.  1  and  2).  That  is,  the  average  star- 
burst  spectral  energy  distribution  is  bluer  at  2rest  =  5-30  pm  than 
those  ULIRGs  with  small  r^/rgg.  Assuming  that  the  average 
spectrum  is  representative  of  the  spectral  energy  distribution  of 
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TABLE  5 
Ice  Absorption 


Object 

6.0  fim 

7*6.0 

3.1  fim 

A.l 

(1) 

(2) 

(3) 

(4) 

(5) 

IRAS  00188-0856 . 

o 

0.9 

O 

1.8 

IRAS  00482-2721 . 

X 

IRAS  03250+1606 . 

X 

0 

0.8 

IRAS  04103-2838 . 

X 

IRAS  08572+3915 . 

X 

X 

IRAS  09039+0503 . 

0 

0.6 

O 

0.8 

IRAS  09116+0334 . 

X 

O 

0.4 

IRAS  09539+0857 . 

X 

X 

IRAS  10378+1108 . 

O 

0.8 

O 

0.6 

IRAS  10485-1447 . 

O 

0.5 

0 

0.8 

IRAS  10494+4424 . 

O 

0.4 

0 

1.0 

IRAS  11095-0238 . 

O 

0.5 

X 

IRAS  11130-2659 . 

O 

0.5 

IRAS  12112+0305 . 

A 

X 

IRAS  12127-1412 . 

O 

0.4 

0 

0.4 

IRAS  12359-0725 . 

X 

0 

0.5 

IRAS  13335-2612 . 

X 

IRAS  14252-1550 . 

X 

0 

0.7 

IRAS  14348-1447 . 

0 

0.9 

o 

0.5 

Arp  220 . 

O 

0.6 

IRAS  16090-0139 . 

O 

0.5 

0 

0.8 

IRAS  16468+5200 . 

A 

X 

IRAS  16487+5447 . 

A 

0 

1.2 

IRAS  17028+5817 . 

A 

0 

0.6 

IRAS  17044+6720 . 

X 

X 

IRAS  21329-2346 . 

A 

0 

1.0 

IRAS  23234+0946 . 

X 

X 

IRAS  23327+2913 . 

X 

X 

IRAS  00091-0738 . 

O 

0.5 

IRAS  00456-2904 . 

X 

IRAS  01004-2237 . 

X 

IRAS  01166-0844 . 

X 

IRAS  01298-0744 . 

O 

0.5 

IRAS  01569-2939 . 

X 

IRAS  02411+0353 . 

X 

IRAS  10190+1322E . 

X 

X 

IRAS  10190+1322W . 

X 

X 

IRAS  11387+4116 . 

X 

X 

IRAS  11506+1331 . 

0 

0.5 

O 

1.0 

IRAS  13509+0442 . 

X 

X 

IRAS  13539+2920 . 

X 

X 

IRAS  14060+2919 . 

X 

X 

IRAS  15206+3342 . 

X 

X 

IRAS  15225+2350 . 

O 

0.5 

O 

0.4 

IRAS  16474+3430 . 

X 

O 

0.8 

IRAS  20414-1651 . 

X 

X 

IRAS  21208-0519 . 

X 

X 

IRAS  22206-2715 . 

X 

IRAS  22491-1808 . 

X 

Notes. — Col.  (1):  Object  name.  Col.  (2):  Detection  or  nondetection  of  the 
6.0  /im  ice  absorption  feature:  (O)  clearly  detected;  (A)  possibly  detected;  (X) 
nondetected.  The  nondetection  includes  two  cases:  (1)  the  absorption  feature  is 
indeed  weak  or  absent  in  high  signal-to-noise  ratios  spectra  or  (2)  the  presence 
of  ice  absorption  is  not  clear  due  to  relatively  large  scatter  in  spectra  or  to  the 
difficulty  in  separating  the  ice  absorption  and  nearby  PAH  emission  features. 
Col.  (3):  Optical  depth  of  the  6.0  fim  ice  absorption  feature.  Ellipses  signify  the 
depth  is  not  measurable.  Since  data  scatter  in  spectra  is  not  the  same  for  individual 
ULIRGs,  there  is  no  fixed  detection  threshold.  The  uncertainty  is  ~0. 1  for  de¬ 
tected  sources.  Col.  (4):  Detection  or  nondetection  of  the  3.1  /mi  ice  absorption 
feature  in  the  2.8-4. 1  /im  spectra  of  Imanishi  et  al.  (2006a):  (O)  clearly  detected; 
(X)  nondetected.  Ellipses  indicate  no  data.  Col.  (5):  Optical  depth  of  the  3.1  /im  ice 
absorption  feature.  Ellipses  signify  the  depth  is  not  measurable  because  (1)  it  is 
weak,  (2)  data  scatter  is  large  in  spectra,  and  (3)  2.8-4. 1  /im  spectra  are  lacking. 
Again  there  is  no  fixed  detection  threshold. 


the  modestly  obscured  starbursts  in  these  ULIRGs,  this  possi¬ 
ble  effect  increases  (not  decreases)  the  observed  Tj'8/r9  7  ratios, 
compared  to  an  uncontaminated  spectrum  of  the  highly  obscured 
core.  Thus,  the  observed  small  t18/79.  7  ratios  in  these  ULIRGs 
are  still  best  explained  by  strong  dust  temperature  gradients 
around  buried  AGNs,  with  centrally  concentrated  energy  source 
geometry. 

Finally,  of  course,  silicate  dust  grains  in  the  cores  of  ULIRGs 
may  have  intrinsically  larger  t[%!t 9  7  ratios  than  the  conserva¬ 
tively  assumed  one.  In  this  case,  if  a  ULIRG  shows  a  decreased 
riVr9  7  rat'°  compared  to  the  intrinsic  one,  due  to  a  dust  tem¬ 
perature  gradient,  but  still  has  an  observed  ratio  similar  to  the 
assumed  one,  then  such  a  ULIRG  would  not  be  classified  as  a 
buried  AGN. 

5.2.4.  Comparison  of  the  Strength  of  Silicate  Absorption 
and  PAH  Emission  Features 

Figure  10a  compares  the  optical  depth  of  the  9.7  pm  silicate 
dust  absorption  (r9  7)  and  the  rest-frame  equivalent  widths  of 
the  6.2  pm  PAH  emission  feature  (EW6.2pah)-  With  the  excep¬ 
tion  of  the  two  outliers  with  both  small  EW6.2pah  and  small  r9  7 
(IRAS  01004-2237  and  17044+6720,  both  weakly  obscured 
AGNs),  there  is  a  weak  anticorrelation,  and  ULIRGs  with  t9.7  ->  ^ 
(i.e.,  buried  AGN  candidates)  tend  to  show  small  EW6.2pah  com¬ 
pared  to  starburst-like  ULIRGs  with  small  r9  7  (<2). 

A  highly  obscured  (Ay  »  30  mag)  normal  starburst  galaxy 
with  mixed  dust/source  geometry  (scenario  A)  predicts  a  rela¬ 
tively  flat  distribution  (i.e.,  relatively  constant  EW6.2pah  values 
with  varying  r9  7),  because  the  PAH  equivalent  width  is  unchanged 
by  dust  extinction  (§  1).  On  the  other  hand,  scenario  B,  where 
(1)  the  PAH  emission  comes  from  modestly  obscured  (Ay  <  20- 
30  mag)  starbursts  exterior  to  the  deeply  obscured  ULIRG  cores 
and  (2)  the  large  r9  7  values  reflect  dust  column  density  toward  the 
buried  AGNs  in  the  cores,  predicts  an  anticorrelation  between 
r9  7  and  EW6.2pah-  As  the  contribution  from  the  AGN-produced 
PAH-free  continuum  to  the  observed  flux  increases,  the  EW6.2pah 
value  decreases.  At  the  same  time,  the  increasing  contribution 
from  the  buried  AGN  emission  (with  centrally  concentrated  en¬ 
ergy  source  geometry)  results  in  the  increase  of  the  observed  r9  7 
value,  compared  to  a  pure  starburst  with  mixed  source/dust  geom¬ 
etry.  The  trend  in  Figure  10a  (see  also  Spoon  et  al.  2007)  sup¬ 
ports  the  scenario  B  for  ULIRGs  with  r9  7  >  2. 

Figure  10b  compares  r9  7  and  the  rest-frame  equivalent  widths 
ofthe  11.3  /im  PAH  emission  feature  (EWiupah)-  Compared  to 
Figure  10a,  the  decreasing  trend  of  EWhjpah  atr97  >  2  appears 
slight  at  best.  In  scenario  B,  the  flux  depression  of  the  AGN- 
produced  continuum  flux  at  Arest  =11.3  pm  is  particularly  large 
due  to  the  strong  9.7  pm  silicate  absorption.  However,  the  1 1 .3  pm 
PAH  emission  flux  is  less  affected  by  this  absorption  because  the 
PAH  emission  comes  from  starburst  regions  outside  the  core. 
Consequently,  the  decrease  of  the  E  Wi  i  jpah  value  with  increas¬ 
ing  t9  7  value  is  small  compared  to  EW6.2pah>  explaining  the 
trends  in  both  Figs.  10a  and  10b. 

Figure  1 1  compares  the  r9  7  values  with  the  1 1 .3-6.2  pm  PAH 
luminosity  ratio.  In  scenario  A,  the  1 1 .3-6.2  pm  flux  ratio  should 
decrease  with  increasing  r9  7,  as  actually  seen  in  normal  starburst 
galaxies  (Brandi  et  al.  2006),  because  the  flux  depression  of  the 
11.3  pm  PAH  emission  inside  the  strong  9.7  pm  absorption  is 
larger.  On  the  other  hand,  in  scenario  B,  no  correlation  is  expected 
because  both  the  6.2  and  11.3  pm  PAH  emission  come  from 
the  modestly  obscured  starbursts,  unrelated  to  the  observed  r9  7 
values.  Figure  11  shows  no  significant  correlation,  supporting 
scenario  B. 


92 


IMANISH1  ET  AL. 


Vol.  171 


Fig.  8. — ( a )  Comparison  of  the  rest-frame  equivalent  width  of  the  6.2  pm  PAH  emission  ( abscissa )  and  6.2  pn i  PAH-to-infrared  luminosity  ratio  {ordinate).  Filled 
triangle :  ULIRGs  classified  optically  as  LINERs.  Open  circle :  ULIRGs  classified  optically  as  H  n  regions.  The  horizontal  dotted  line  indicates  the  ratio  for  normal 
starburst  galaxies  with  modest  dust  obscuration  (T6.2PAh/^ir  —  3.4  x  10-3;  see  §  5. 1).  In  the  upper  left  small  area  separated  by  the  dashed  line,  there  are  no  data  points. 
(b)  The  abscissa  is  the  rest-frame  equivalent  width  of  the  11.3  pm  PAH  emission,  and  the  ordinate  is  the  11.3  pm  PAH-to-infrared  luminosity  ratio.  The  horizontal  dotted 
line  indicates  the  ratio  for  normal  starburst  galaxies  with  modest  dust  obscuration  (£ii.3Pah/^ir  =  1  -4  x  10-3;  see  §5.1).  The  larger  scatter  of  L\ i.3pah/^ir  in  this  figure 
than  for  L6.2pah/Tir  in  (a)  may  be  due  to  intrinsically  larger  scatter  of  the  11.3  pm  PAH  strengths  in  normal  starburst  galaxies  (see  §  5.2.1). 


5.2.5.  Is  the  Buried  AGN  Scenario  Unique? 

ULIRGs  with  large  r'g  7  (>2)  and  small  t[ 8/t9'  7  ratios  are  nat¬ 
urally  explained  by  the  presence  of  buried  AGNs  (scenario  B  ). 
We  here  briefly  comment  on  two  alternative  scenarios  that  are 
possible  but  we  believe  are  unlikely. 

First,  in  a  normal  starburst  with  mixed  dust/source  geometry, 
if  each  star  is  deeply  buried  in  a  large  amount  of  circumstellar 
dust  cocoon,  such  as  a  protostar,  then  a  strong  dust  absorption 
feature  (t^  7  >  2)  could  be  produced.  However,  the  timescale 
of  such  a  protostar  phase  is  only  <106  yr  (Lada  1987),  which 
is  much  less  than  the  typical  timescale  of  a  normal  starburst 
(~  1 07- 1 08  yr;  Moorwood  1 996)  or  even  the  maximum  starburst 
(~  1 07- 1 0s  yr;  Elmegreen  1999),  or  a  ULIRG  timescale  (~108- 
109  yr;  Murphy  et  al.  1996).  It  is  therefore  very  unlikely  that  all 
dust-heating  stars  are  protostars  in  a  significant  fraction  of  ULIRGs. 
Furthermore,  this  scenario  produces  no  strong  dust  temperature 
gradient  and  so  cannot  reduce  the  observed  r[&/Tg  7  ratio. 

Second,  it  may  be  argued  that  exceptionally  centrally  concen¬ 
trated  starbursts  remain  a  possibility  (Fig.  le).  However,  while 
starburst  galaxies  with  Zir  =  109,7-1011,6  L0  are  well  explained 
by  the  mixed  dust/source  geometry  ( Brandi  et  al.  2006;  see  §  5.2.4 
of  our  paper),  there  is  no  reason  why  starbursts  at  the  cores  of 
ULIRGs  (Z|r  >  1012  Lq)  suddenly  change  their  properties  to 
exceptionally  centrally  concentrated  nature  (Fig.  le).  Further¬ 
more,  Soifer  et  al.  (2000)  studied  the  cores  of  seven  nearby 
ULIRGs  and  found  that  the  observed  emission  surface  bright¬ 
nesses  in  the  majority  (six  out  of  seven)  of  these  ULIRGs’  cores 
are  very  high  (^1013  L0  kpc~2).  Even  if  the  energy  sources  are 
uniformly  distributed  over  the  cores  (<500  pc  in  size),  the  ob¬ 
served  surface  brightnesses  are  close  to  the  maximum  values 
seen  in  the  cores  of  Galactic  H  n  regions  (Soifer  et  al.  2000).  If  an 
exceptionally  centrally  concentrated  starburst  were  the  luminosity 
source  in  the  core  of  a  ULIRG,  the  emission  surface  brightness  of 
the  starburst  would  have  to  be  extremely  high,  >1014  L0  kpc~2. 
Among  known  star-forming  regions,  such  a  high  surface  bright¬ 
ness  could  only  barely  be  achieved  with  a  super  star  cluster,  where 
many  young  stars  are  formed  within  an  area  of  less  than  a  few 
parsec  (Goijian  et  al.  2001).  The  absolute  luminosity  of  such  super 
star  clusters  is,  however,  only  of  order  109— 1010  L@,  and  so  the 


super  star  cluster  scenario  may  work  for  the  centrally  concen¬ 
trated  energy  sources  of  low-luminosity  (LiR  ~  1010  L0)  galaxy 
nuclei  (Imanishi  &  Nakanishi  2006).  For  the  centrally  concen¬ 
trated  energy  sources  at  ULIRG  cores,  both  high  emission  surface 
brightness  and  large  absolute  luminosity  (^1012  Lq)  are  required, 
making  the  buried  AGN  scenario  the  most  plausible  choice  (Soifer 
et  al.  2003,  their  §  6).  In  fact,  known  AGN  populations  (i.e.,  qua¬ 
sars)  meet  these  requirements  without  difficulty,  while  for  the  star- 
burst  scenario,  of  order  100-1000  super  star  clusters  must  be 
concentrated  in  a  very  compact  region  and  such  a  phenomenon 
has  not  been  clearly  confirmed  from  observations.  We  thus  believe 
that  buried  AGNs  are  the  most  natural  explanation  for  the  cen¬ 
trally  concentrated  energy  sources  at  ULIRGs’  cores.  The  pres¬ 
ence  of  luminous  buried  AGNs  in  many  nearby  ULIRGs  agrees 
with  the  suggestion  by  Takeuchi  et  al.  (2003)  that  the  observed 
lump  of  the  local  infrared  60  pm  luminosity  function  at  the  highest 
luminosity  end  is  due  to  the  AGN  contribution. 

5.2.6.  Combination  of  Energy  Diagnostic  Methods 

Table  6  (col.  [7])  summarizes  the  strengths  of  the  detected 
buried  AGN  signatures  in  Spitzer  IRS  5-35  pm  spectra  based 
on  the  three  methods:  (1)  low  PAH  equivalent  width,  (2)  large 
Tgj  value,  and  (3)  small  t{8 /tj  7  and  t^/t^  7  ratio.  When  buried 
AGN  signatures  in  individual  ULIRGs  are  consistently  found 
using  all  or  most  of  the  methods,  the  ULIRGs  are  classified  as 
very  strong  buried  AGN  candidates,  marked  with  open  double 
circles.  When  buried  AGNs  signatures  are  seen  in  the  first  method 
(  particularly  based  on  the  observed  EW6.2pah  value)  and/or  in  the 
third  method,  then  the  ULIRGs  are  classified  as  strong  AGN  can¬ 
didates  ( open  circles).  When  the  signatures  are  detected  only  in 
the  second  method  and/or  low  observed  EWhqpah  or  corrected 
EW6.2pah  values  (first  method),  the  ULIRGs  are  classified  as  pos¬ 
sible  buried  AGN  candidates  ( open  triangles ),  because  (1)  a  nor¬ 
mal  starburst  nucleus  obscured  by  foreground  dust  in  the  host 
galaxy  (Fig.  Id)  cannot  be  ruled  out  completely  and  (2)  the  11.3  pm 
PAH  emission  strengths  in  normal  starbursts  have  intrinsically 
larger  scatter  (§  5.2.1).  The  fraction  of  ULIRGs  with  strong  bur¬ 
ied  AGN  signatures  is  16/48  (33%;  LINER  +  H  n  region),  10/28 
(36%;  LINER),  and  6/20  (30%;  H  n  regions).  When  ULIRGs  with 
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Fig.  9. — {a)  Comparison  of  the  rest-frame  equivalent  widths  of  PAH  emission  features  in  nm.  The  abscissa  is  the  6.2  pm  PAH,  and  the  ordinate  is  the  7.7  pn l  PAH 
emission.  Filled  triangle :  ULIRGs  classified  optically  as  LINERs.  Open  circle :  ULIRGs  classified  optically  as  H  n  regions.  In  the  upper  left  small  area  separated  with  the 
dashed  line,  there  are  no  data  points.  ( b )  The  abscissa  is  the  6.2  pm  PAH,  and  the  ordinate  is  the  11.3  pm  PAH  emission,  (c)  The  abscissa  is  the  7.7  pm  PAH,  and  the  ordi¬ 
nate  is  the  1 1 .3  pm  PAH  emission.  ( d )  The  abscissa  is  the  3.3  pm  PAH,  and  the  ordinate  is  the  6.2  pm  PAH  emission,  (e)  The  abscissa  is  the  3.3  pm  PAH,  and  the  ordinate 
is  the  7.7  pm  PAH  emission.  (/)  The  abscissa  is  the  3.3  pm  PAH,  and  the  ordinate  is  the  1 1.3  pm  PAH  emission.  In  the  last  three  plots,  the  average  EW3 jpah  value  of 
both  nuclei  is  plotted  for  IRAS  12112+0305  and  14348—1447  (see  Table  3). 


possible  buried  AGN  signatures  are  included,  the  fraction  in¬ 
creases  to  26/48  (54%;  LINER  +  H  n  region),  19/28  (68%; 
LINER),  and  7/20  (35%;  H  n  region).  Our  careful  look  at  the 
infrared  5-35  pm  spectra  suggests  that  luminous  buried  AGNs  re¬ 
side  in  a  much  higher  fraction  of  optical  non-Seyfert  ULIRGs  than 
previous  estimates  (Taniguchi  et  al.  1999;  Lutz  et  al.  1999). 

We  stress  that  the  above  fraction  of  luminous  buried  AGNs  is 
only  for  ULIRGs  classified  optically  as  LINERs  and  H  n  regions 
(i.e.,  non-Seyferts),  which  comprise  70%  of  ULIRGs.  Given  that 
30%  of  ULIRGs  show  optical  Seyfert  signatures,  indicative  of 
luminous  AGNs  surrounded  by  dusty  tori,  and  that  the  AGNs  are 
estimated  to  be  energetically  dominant  in  these  optical  Seyfert 
ULIRGs  (Veilleux  et  al.  1999b),  the  total  fraction  of  ULIRGs 
possessing  luminous  AGNs  is  at  least  50%,  and  possibly  >65%, 


in  the  whole  UL1RG  sample  at  z  <  0.15,  including  both  optical 
Seyfert  and  non-Seyfert  objects. 

Finally,  in  all  of  the  methods  employed  by  us  (§  5.2.1 -5.2.3 
and  5.2.6),  ULIRGs  classified  optically  as  LINERs  tend  to  dis¬ 
play  signatures  of  luminous  buried  AGNs  more  frequently  than 
those  classified  optically  as  H  n  regions.  For  starburst/buried 
AGN  composite  ULIRGs  (Fig.  lc),  it  is  likely  that  the  optical 
LINER  or  H  n  region  classifications  are  largely  affected  by  the 
properties  of  the  modestly  obscured  starbursts  at  the  exteriors  of 
the  ULIRG  cores,  rather  than  buried  AGN-related  emission,  be¬ 
cause  optical  observations  can  probe  only  the  surfaces  of  dusty 
objects.  If  shocks  or  superwinds  are  important  in  the  modestly 
obscured  starbursts,  the  optical  classification  will  be  LINERs.  If 
emission  from  H  n  regions  is  dominant,  the  ULIRGs  will  be 


TABLE  6 

Buried  AGN  Signatures 


Object 

(i) 

Spitzer  IRS  5- 

-35  fim 

3-4  fim 
(8) 

EW62pah 

(2) 

EW7  7pah 
(3) 

EW11  3PAH 

(4) 

r9.7 

(5) 

T-Gradient 

(6) 

Total 

(7) 

IRAS  00188-0856 . 

.  O 

X 

X 

0 

0 

© 

0 

IRAS  00482-2721 . 

.  X 

X 

X 

0 

X 

A 

IRAS  03250+1606 . 

.  X 

X 

X 

X 

X 

X 

0 

IRAS  04103-2838 . 

.  O 

X 

0 

X 

X 

O 

IRAS  08572+3915 . 

.  O 

0 

0 

0 

0 

© 

© 

IRAS  09039+0503 . 

.  X 

X 

X 

0 

X 

A 

0 

IRAS  09116+0334 . 

.  X 

X 

X 

X 

X 

X 

0 

IRAS  09539+0857 . 

.  0 

X 

X 

0 

X 

O 

X 

IRAS  10378+1108 . 

.  O 

X 

X 

0 

O 

0 

0 

IRAS  10485-1447 . 

.  X  (O) 

X 

0 

0 

A 

A 

0 

IRAS  10494+4424 . 

.  X 

X 

X 

X 

X 

X 

0 

IRAS  11095-0238 . 

.  O 

X 

0 

0 

0 

© 

X 

IRAS  11130-2659 . 

.  O 

X 

X 

0 

X 

0 

X 

IRAS  12112+0305 . 

.  X 

X 

X 

X 

X 

X 

X 

IRAS  12127-1412 . 

.  O 

0 

0 

0 

0 

© 

© 

IRAS  12359-0725 . 

.  X 

X 

0 

X 

A 

A 

0 

IRAS  13335-2612 . 

.  X 

X 

X 

X 

X 

X 

IRAS  14252-1550 . 

.  X 

X 

X 

X 

X 

X 

0 

IRAS  14348-1447 . 

.  X  (O) 

X 

X 

X 

X 

A 

0 

Arp  220 . 

.  x  (O) 

X 

X 

0 

X 

A 

X 

IRAS  16090-0139 . 

.  O 

X 

X 

0 

A 

O 

0 

IRAS  16468+5200 . 

.  X 

X 

X 

0 

A 

A 

X 

IRAS  16487+5447 . 

.  X 

X 

X 

X 

X 

X 

0 

IRAS  17028+5817 . 

.  X 

X 

X 

X 

X 

X 

0 

IRAS  17044+6720 . 

.  O 

0 

0 

X  (.  .  ,)a 

O 

© 

© 

IRAS  21329-2346 . 

.  X 

X 

X 

0 

X 

A 

0 

IRAS  23234+0946 . 

.  X 

X 

X 

X 

X 

X 

X 

IRAS  23327+2913 . 

.  X 

X 

0 

X 

X 

A 

X 

IRAS  00091-0738 . 

.  0 

X 

X 

0 

X 

O 

IRAS  00456-2904 . 

.  X 

X 

X 

X 

X 

X 

IRAS  01004-2237 . 

.  0 

0 

0 

X  (. .  .y 

X 

© 

IRAS  01166-0844 . 

.  0 

X 

0 

0 

0 

© 

IRAS  01298-0744 . 

.  0 

X 

0 

0 

X 

0 

IRAS  01569-2939 . 

.  0 

X 

X 

O 

X 

O 

IRAS  02411+0353 . 

.  X 

X 

X 

X 

X 

X 

IRAS  10190+1322 . 

.  xb 

xb 

xb 

X 

X 

X 

X 

IRAS  11387+4116 . 

.  X 

X 

X 

X 

X 

X 

X 

IRAS  11506+1331  . 

.  X  (O) 

X 

X 

O 

X 

A 

O 

IRAS  13509+0442 . 

.  X 

X 

X 

X 

X 

X 

X 

IRAS  13539+2920 . 

.  X 

X 

X 

X 

X 

X 

X 

IRAS  14060+2919 . 

.  X 

X 

X 

X 

X 

X 

X 

IRAS  15206+3342 . 

.  X 

X 

X 

X 

X 

X 

X 

IRAS  15225+2350 . 

.  0 

X 

0 

O 

X 

0 

O 

IRAS  16474+3430 . 

.  X 

X 

X 

X 

X 

X 

O 

IRAS  20414-1651 . 

.  X 

X 

X 

X 

X 

X 

X 

IRAS  21208-0519 . 

.  X 

X 

X 

X 

X 

X 

X 

IRAS  22206-2715 . 

.  X 

X 

X 

X 

X 

X 

IRAS  22491-1808 . 

.  X 

X 

X 

X 

X 

X 

Notes. — Col.  (1):  Object  name.  Col.  (2):  Buried  AGN  signatures  based  on  the  low  equivalent  width  of  the  6. 2  pm  PAH  emission  (EW6.2pah  <  180  nm;  §  5.2.1):  (O) 
present;  (X)  none.  When  the  corrected  EW^pah  value  becomes  less  than  180  nm,  the  symbol  “O”  is  shown  in  parentheses.  Col.  (3):  Buried  AGN  signatures  based  on 
the  low  equivalent  width  of  the  7.7  pm  PAH  emission  (EW7.7PAH  <  230  nm;  §  5.2.1):  (O)  present;  (X)  none;  Col.  (4):  Buried  AGN  signatures  based  on  the  low 
equivalent  width  of  the  1 1 .3  pm  PAH  emission  (EWhjpah  <  200  nm;  §  5.2.1):  (O)  present;  (X)  none.  Col.  (5):  Buried  AGN  signatures  based  on  the  large  Tq  7  value 
(>2;  §  5.2.2):  (O)  present;  (X)  none.  Col.  (6):  Buried  AGN  signatures  based  on  the  small  7  and  Tig/79 7  ratios:  (O)  present;  (A)  possibly  present;  (X)  none. 
Considering  the  uncertainty  of  the  ratio  with  <10%  (Table  4),  ULIRGs  showing  7  <  0.27  and  7  <  0.32  are  taken  as  sources  with  detectable  signatures  of 
strong  dust  temperature  gradients  (§  5.2.3).  Col.  (7):  Buried  AGN  signatures  from  combined  methods  in  Col.  (2)-(6)  based  on  Spitzer  IRS  spectra.  ©:  very  strong:  (O) 
strong;  (A)  possibly  present;  (X)  none.  Please  see  §  5.2.6  for  more  details.  Col.  (7):  Buried  AGN  signatures  based  on  infrared  3-4  pm  spectroscopic  diagnostic  by 
Imanishi  et  al.  (2006a):  (©)  very  strong;  (O)  strong;  (X)  none;  (. . .)  no  data. 

a  The  small  7-9  7  value  is  due  to  weak  obscuration  of  a  luminous  AGN,  rather  than  the  normal  starburst  nature  with  mixed  dust /source  geometry. 

PAH  equivalent  widths  are  derived  for  both  nuclei  separately.  No  AGN  signatures  are  found  in  both  nuclei. 
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Fig.  10. — (a)  Comparison  of  the  strength  of  the  9.7  fim  silicate  dust  absorption  (7-9  7)  ( abscissa )  and  the  rest-frame  equivalent  widths  of  the  6.2  fim  PAH  emission 
feature  (ordinate).  Filled  triangle :  ULIRGs  classified  optically  as  LINERs.  Open  circle :  ULIRGs  classified  optically  as  H  11  regions.  In  the  upper  right  small  area 
separated  with  the  dashed  line,  there  are  no  data  points,  (b)  The  ordinate  is  the  rest- frame  equivalent  width  of  the  11.3  fin 1  PAH  emission  feature.  In  both  plots,  the  two 
outliers  IRAS  01004—2237  and  17044+6720  (weakly  obscured  AGNs;  see  §  5.2.4)  are  indicated  as  “IR01004”  and  “IR17044,”  respectively.  For  IRAS  10190+1322, 
the  average  value  of  both  nuclei  (see  Table  3)  are  plotted  for  PAH  equivalent  widths  in  the  ordinate. 


classed  as  H  n  regions.  In  this  case,  it  is  not  obvious  why  two 
phenomena,  a  luminous  buried  AGN  at  the  center  and  shocks 
in  the  surface  starburst,  are  correlated.  Veilleux  et  al.  (1995, 
1999a)  found  that  the  emission  probed  in  the  optical  is  dustier  in 
LINER  ULIRGs  than  in  H  n  region  ULIRGs.  In  a  dusty  star- 
burst,  shock-related  emission  can  be  relatively  important  in  the 
optical  compared  to  the  emission  from  H  n  regions  themselves, 
resulting  in  optical  LINER  classification.  When  a  luminous  AGN 
is  placed  at  the  center  of  a  less  dusty  starburst,  classified  optically 
as  an  H  n  region,  the  AGN  emission  is  more  easily  detectable  in 
the  optical  spectrum,  making  such  an  object  an  optical  Seyfert.  On 
the  other  hand,  when  a  luminous  AGN  is  placed  at  the  center  of  a 
dusty  starburst  classified  optically  as  a  LINER,  the  AGN  emission 


Fig.  11. — Comparison  of  the  strength  of  the  9.7  fin i  silicate  dust  absorption 
(79  7)  (abscissa)  and  11.3-6.2  fin 1  PAH  luminosity  ratio  (ordinate).  Filled  tri¬ 
angle :  ULIRGs  classified  optically  as  LINERs.  Open  circle :  ULIRGs  classified 
optically  as  H  11  regions.  In  the  upper  left  small  area  separated  with  the  dashed 
line,  there  are  no  data  points.  For  IRAS  10190+1322,  the  11.3-6.2  fim  PAH 
luminosity  ratio  in  the  ordinate  is  derived  by  combining  spectra  of  both  nuclei. 
IRAS  08572+3915  is  not  plotted,  because  both  the  6.2  and  1 1.3  fin 1  PAH  emis¬ 
sion  luminosities  are  upper  limits  (see  Table  3). 


is  more  elusive  in  the  optical  spectrum,  so  that  such  an  object 
is  classified  as  an  optical  non-Seyfert.  The  observed  higher  frac¬ 
tion  of  optically  elusive  buried  AGN  s  in  optical  LINER  ULIRGs 
compared  to  H  n  region  ULIRGs  can  qualitatively  be  explained 
by  this  scenario. 

5.2.7.  Absorption-corrected  Luminosities  of  Buried  AGNs 

For  ULIRGs  that  show  buried  AGN  signatures  and  also  small 
PAH  equivalent  widths,  the  observed  fluxes  are  largely  ascribed 
to  AGN-heated  dust  continuum  emission.  The  resulting  luminos¬ 
ity  is  conserved  at  each  temperature  (Fig.  2).  We  can  quantitatively 
estimate  the  absorption-corrected  intrinsic  dust  luminosity  heated 
by  the  AGN  (whole  infrared),  or  the  absorption-corrected  AGN 
energetic  radiation  luminosity  (X-ray-UV-optical),  based  on  the 
observed  fluxes  at  Arest  ~  10  pm  and  the  dust  column  density  to¬ 
ward  the  10  pm  continuum-emitting  regions  inferred  from  t, 97. 
Since  in  a  buried  AGN,  the  10  /(in  continuum  emission  regions 
are  dominated  by  fixed  layers  (§  5.2.3;  Fig.  2),  the  foreground 
screen  dust  absorption  model  is  applicable,  and  the  dust  extinc¬ 
tion  corrections  are  straightforward. 

To  do  this,  we  assume  that  r9  7  and  the  extinction  at  2rest  =  8 
or  13  pm  continuum  just  outside  the  9.7  pm  silicate  feature 
G+ont)  are  related  with  Tg  7/Acont  ~  2.3  (Rieke  &  Lebofsky  1985). 
Table  7  summarizes  the  absorption-corrected  AGN  luminosities 
for  selected  ULIRGs  with  both  strong  buried  AGN  signatures 
and  low  PAH  equivalent  widths.  The  flux  attenuation  of  the  8  or 
13  pm  continuum  outside  the  9.7  pm  silicate  feature  ranges 
from  a  factor  of  1 .2  (IRAS  01004-2237;  t£7  ~  0.4)  to  5  (IRAS 
01298—0744;  r9  7  ~  4).  The  absorption-corrected  AGN  lumi¬ 
nosities  are  (l-5.5)xl045  ergs  s~ 1 ,  which  are  of  order  ~  1 0 1 2  LQ 
(=3.8  x  1045  ergs  s_1). 

The  ambiguity  of  the  derived  absorption-corrected  AGN 
luminosities  come  from  (1)  absolute  flux  calibration  error  of  the 
Spitzer  IRS  spectra  and  slit  loss  of  the  compact  AGN  emission; 
(2)  possibly  different  dust  extinction  curves  from  the  Galactic 
one  that  we  assume;  and  (3)  a  small  contribution  from  PAH- 
emitting,  modestly  obscured  starbursts,  which  decreases  the  r9'  7 
values  and  increases  the  continuum  flux  level,  compared  to  pure 
emission  from  the  buried  AGNs.  The  first  error  is  not  significant, 
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TABLE  7 

Absorption-corrected  AGN  Luminosity 


Object 

L(AGN) 
(1045  ergs  s-1) 

U) 

(2) 

IRAS 

IRAS 

IRAS 

IRAS 

IRAS 

IRAS 

IRAS 

IRAS 

IRAS 

IRAS 

IRAS 

IRAS 

IRAS 

IRAS 

IRAS 


00188-0856 .  1.5 

04103-2838 .  1.5 

08572+3915 .  5.5 

10378+1108 .  1.0 

10485-1447 .  1.0 

11095-0238 .  2.5 

11130-2659 .  1.0 

12127-1412 .  2.5 

16090-0139 .  2.0 

17044+6720 .  2.0 

00091-0738 .  1.5 

01004-2237 .  1.5 

01166-0844 .  1.0 

01298-0744 .  3.5 

15225+2350 .  1.5 


Notes. — Col.  (1):  Object  name.  Col.  (2):  Absorption- 
corrected  intrinsic  luminosity  of  the  AGN  energetic  radiation 
( X-ray— UV-optical)  in  units  of  1045  ergs  s-1. 


approximately  20%  or  so.  The  second  uncertainty  is  more  diffi¬ 
cult  to  quantify,  but  our  adoption  of  79  7  reduces  the  effects  of  dust 
grain  size  change  (§  5.2.2).  For  the  third  source  of  uncertainty, 
the  effect  of  the  decreased  r9  7  values  is  more  important  than  that 
of  the  increased  continuum  flux  level,  because  the  AGN-produced 
flux  at  Arest  ~  9.7  jun  is  more  highly  attenuated.  Assume  for  ex¬ 
ample  that  the  modestly  obscured  starburst  emission  contributes 
1 0%  of  the  observed  flux  at  2rest  =  8  and  1 3  pm,  outside  the  9.7  /im 
silicate  feature.  The  buried  AGN-produced  continuum  flux  is 
then  only  10%  smaller  than  the  observed  flux.  For  ULIRGs  with 
Tgj  >  2  (i.e.,  buried  AGN  candidates  we  are  now  discussing), 
since  the  attenuation  of  the  observed  flux  at  2rest  =  9.7  pm  is  a 
factor  of  >7  compared  to  the  nearby  continuum,  more  than  50% 
of  the  continuum  at  2rest  =  9.7  71m  is  produced  by  the  starburst 
emission.  The  true,  uncontaminated  t9  7  value  toward  the  buried 
AGN  is  larger  by  a  factor  of  >0.7  than  the  observed  value,  so  that 
the  correction  factor  of  dust  extinction  for  the  AGN -heated  dust 
continuum  is  larger  by  >30%  (r9'  7A4co nt  ~  2.3).  Therefore,  the 
small  starburst  contamination  generally  results  in  the  under¬ 
estimate  of  the  intrinsic  absorption-corrected  AGN  luminosity. 
Despite  these  uncertainties,  the  rough  agreement  between  the 
absorption-corrected  buried  AGN  luminosities  and  the  observed 
infrared  luminosities  of  these  ULIRGs  [(1-3)  x  1012  Z0]  makes  it 
reasonable  for  us  to  argue  that  the  putative  AGNs  could  quanti¬ 
tatively  account  for  the  bulk  of  the  observed  infrared  luminosities, 
at  least  for  ULIRGs  with  strong  buried  AGN  signatures  in  Spitzer 
IRS  spectra. 

5.3.  Comparison  with  Infrared  3-4  pm  Spectra 

5.3.1.  Ice  Absorption  Feature 

H20  ice  absorption  is  also  present  at  3. 1  pm  (due  to  a  stretch¬ 
ing  mode)  and  has  been  detected  in  3-4  pm  spectra  of  many 
ULIRGs  in  this  sample  (Imanishi  &  Maloney  2003;  Imanishi 
et  al.  2006a;  Risaliti  et  al.  2006).  Comparison  of  the  3. 1  and  6.0  pm 
ice  absorption  features  can  provide  a  better  understanding  of  the 
properties  of  ice  in  ULIRGs. 

Column  (4)  of  Table  5  summarizes  the  detection  or  non¬ 
detection  of  the  3.1  pm  ice  feature.  For  the  detections  its  optical 
depth  is  shown  in  column  (5).  Many  ULIRGs  have  clear  detec¬ 
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Fig.  12. — Comparison  of  the  optical  depths  of  the  H2O  ice  absorption  fea¬ 
tures  at  3. 1  /rm  ( abscissa )  and  6.0  jim  {ordinate).  The  dashed  line  is  the  expected 
ratio  from  a  laboratory  H2O  ice  sample  at  10  K  (Keane  et  al.  2001). 


tions  of  both  the  3.1  and  6.0  pm  ice  absorption  features,  re¬ 
inforcing  the  picture  that  the  energy  sources  in  many  ULIRGs 
are  obscured  by  a  large  amount  of  ice-covered  dust  grains. 

Figure  12  compares  the  optical  depths  of  the  3. 1  pm  (t31)  and 
6.0  pm  (r6  o)  H20  ice  absorption  features,  for  ULIRGs  clearly 
showing  both  features.  The  T6.o/t3.i  ratios  are  several  times  larger 
than  expected  from  laboratory  data.  Systematically  larger  t^.o/t^a 
ratios  compared  to  the  laboratory  prediction  have  previously  been 
found  in  Galactic  obscured  stars,  but  the  reason  is  unclear  (Gibb 
et  al.  2000;  Keane  et  al.  2001).  This  trend  may  be  common  in 
astronomical  obscured  sources.  Our  poor  understanding  of  ice 
absorption  features  hampers  the  search  for  buried  AGNs  with  cen¬ 
trally  concentrated  energy  source  geometry  based  on  small  76.o /t3.  i 
ratios,  as  was  done  based  on  the  small  rlVT9.7  ratio  (§  5.2.3). 

5.3.2.  PAH  Emission  Features 

The  3.3  pm  PAH  emission  feature  assigned  to  the  aromatic 
CH  stretch  (Allamandola  et  al.  1989)  was  detected  in  the  ma¬ 
jority  of  this  ULIRG  sample  (Imanishi  et  al.  2006a).  Figures  9  d, 
9e,  and  9 f  compare  the  rest-frame  equivalent  widths  of  the 

3.3  pm  PAH  feature  (EW3.3pah)  with  those  of  the  6.2,  7.7,  and 

11.3  pm  PAH  features. 

In  Figure  9,  the  scatter  is  generally  large  in  plots  comparing 
PAH  features  at  widely  separated  wavelengths.  Contribution 
from  PAH-free  AGN  continuum  to  an  observed  flux  can  vary  to 
a  larger  degree  with  increasing  wavelength  separation,  explain¬ 
ing  the  trend  in  Figure  9. 

5.3.3.  Support  for  the  Strong  Dust  Temperature  Gradient 

In  a  buried  AGN  with  a  strong  dust  temperature  gradient,  dust 
extinction  toward  the  3-4  pm  continuum  emission  regions,  es¬ 
timated  using  3-4  pm  data,  A  y(3  pm),  should  be  larger  than  that 
toward  10  pm  continuum  emission  region  estimated  from  the 
Tg  7  value  (Fig.  2).  For  ULIRGs  with  small  contamination  from 
PAH-emitting  starbursts  to  their  3-4  pm  spectra  (small  EW3  3pah 
values),  Av (3  pm)  can  be  estimated  with  little  ambiguity,  from 
absorption  optical  depths  at  3.4  pm  (r3  4)  by  bare  carbonaceous 
dust  and  at  3.1  pm  (r3.i)  by  ice-covered  dust  deep  inside  mo¬ 
lecular  clouds. 

IRAS  08572+3915, 1 2 127- 14 1 2,  and  1 7044+6720  show  3- 
4  pm  spectra  with  very  small  PAH  contaminations.  The  r3  4  val¬ 
ues  are  0.8,  0.35,  and  0.15  for  IRAS  08572+3915, 12127-1412, 
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Fig.  13. — Comparison  of  the  strength  of  the  9.7  pm  silicate  dust  absorption  and  infrared  continuum  colors,  (a)  The  abscissa  is  the  t 9  7  value,  and  the  ordinate  is  the 
5.4-35  pm  flux  ratio  measured  from  the  Spitzer  IRS  spectra  in  Figs.  3  and  4.  Filled  triangle :  ULIRGs  classified  optically  as  LINERs.  Open  circle :  ULIRGs  classified 
optically  as  H 11  regions.  In  the  upper  right  small  area  separated  with  the  dashed  line,  there  are  no  data  points,  (b)  The  ordinate  is  the  IRAS  25-60  pm  flux  ratio,  taken  from 
Table  1.  The  horizontal  dotted  line  indicates  the  separation  between  warm  and  cool  ULIRGs  (Sanders  et  al.  1988b). 


and  17044+6720,  respectively  (Imanishi  et  al.  2006a).  The  T3  4 
value  reflects  the  column  density  of  bare  carbonaceous  dust  (with¬ 
out  an  ice  mantle;  Pendleton  et  al.  1994;  Imanishi  et  al.  1996; 
Rawlings  et  al.  2003).  In  the  Galaxy,  the  r9'  7/r3.4  ratio  is  surpris¬ 
ingly  similar  in  different  directions  (7-5  7/73.4  =  14-17;  Roche  & 
Aitken  1984,  1985;  Pendleton  et  al.  1994).  Thus,  a  r9  7/73.4  ratio 
substantially  smaller  than  14-17  can  be  a  signature  of  a  strong 
dust  temperature  gradient.  Adopting  the  r9  7  values  shown  in 
Table  4,  we  obtain  7-57/73  4  =  5,  7,  and  12  for  IRAS  08572+ 
3915,  12127—1412,  and  17044+6720,  respectively.  These  three 
ULIRGs  have  previously  been  classified  as  very  strong  buried 
AGN  candidates,  based  on  Spitzer  IRS  5-35  pm  spectra  (§  5.2.6). 
The  small  r9  7/73.4  ratios  support  the  buried  AGN  hypothesis. 

For  IRAS  08572+3915  and  17044+6720,  no  3.1  pm  ab¬ 
sorption  by  ice-covered  dust  is  detected  (Imanishi  et  al.  2006a), 
suggesting  that  most  of  the  dust  toward  the  3-4  pm  continuum 
emission  sources  is  bare,  without  an  ice  mantle.  However,  for 
IRAS  12127 — 1412,  the  3. 1  pm  absorption  feature  is  found  with 
73.1  ~  0.4,  suggesting  a  large  amount  of  ice-covered  dust,  in  ad¬ 
dition  to  bare  dust  grains.  Since  the  9.7  pm  absorption  is  seen 
whether  the  silicate  dust  is  bare  or  ice-covered,  adding  the  column 
density  of  ice-covered  dust  to  that  of  bare  carbonaceous  dust  grains 
will  raise  Ay(3  pm)  even  further,  strengthening  the  argument  for 
a  strong  dust  temperature  gradient  in  IRAS  12127—1412. 

5.3.4.  Classification:  AGN  or  Starburst? 

Table  6  (col.  [8])  summarizes  the  detection  or  nondetection, 
and  the  detection  significance  in  the  case  of  detection,  of  bur¬ 
ied  AGN  signatures  based  on  ground-based  3-4  pm  spectra 
(Imanishi  et  al.  2006a).  The  trend  of  buried  AGN  detection  is 
generally  consistent  in  individual  ULIRGs,  in  that  (1)  very  strong 
buried  AGN  candidates  at  3-4  pm  are  also  very  strong  buried 
AGN  candidates  in  Spitzer  IRS  5-35  pm  spectra,  and  (2)  ULIRGs 
with  no  AGN  signatures  at  3-4  pm  usually  show  no  AGN  signa¬ 
tures  at  5-35  pm.  There  are  some  ULIRGs  that  have  detectable 
AGN  signatures  at  5-35  pm,  but  not  at  3-4  pm,  which  are  ex¬ 
plained  by  a  larger  contamination  from  modestly  obscured  starburst 
emission  at  a  shorter  wavelength  (Imanishi  et  al.  2004;  Imanishi 
2006). 

Finally,  Imanishi  et  al.  (2006a)  have  previously  argued,  based 
on  3-4  pm  spectra,  that  the  fraction  of  ULIRGs  showing  lu¬ 


minous  buried  AGN  signatures  is  higher  in  LINER  ULIRGs  than 
H  11  region  ULIRGs  in  this  sample.  This  trend  is  exactly  what 
we  have  seen  in  the  Spitzer  IRS  5-35  pm  spectra  discussed  here 
(§  5.2.1-5.2.6). 

5.4.  Buried  AGNs  and  Infrared  Colors 

Figure  13  compares  the  strength  of  the  9.7  pm  silicate  absorp¬ 
tion  feature  (r97)  and  infrared  colors,  Spitzer  IRS  5.4-35  pm 
flux  ratio  (Fig.  13a),  and  IRAS  25-60  pm  flux  ratio  (Fig.  13b). 
The  infrared  colors  are  independent  of  r9  7  values,  with  no  clear 
positive  or  negative  correlation  found. 

The  IRAS  25-60  pm  flux  ratio  is  often  used  to  determine 
whether  ULIRGs  show  cool  (./25//60  <  0.2)  or  warm  (>0.2)  far- 
infrared  colors  (Sanders  et  al.  1988b).  Although  many  classical 
Seyfert-type  AGNs  surrounded  by  dusty  tori  show  warm  far- 
infrared  colors  (de  Grijp  et  al.  1987),  the  bulk  of  ULIRGs  possess¬ 
ing  luminous  buried  AGN  candidates  (r9  7  >  2)  belong  to  cool 
sources  (Fig.  13b).  As  already  argued  by  Imanishi  et  al.  (2006a) 
it  is  quite  reasonable  that  buried  AGNs  tend  to  show  cooler  in¬ 
frared  colors  than  optical  Seyfert  AGNs.  A  larger  amount  of  dust 
around  a  central  AGN  in  a  buried  AGN  compared  to  a  Seyfert 
AGN  surrounded  by  a  torus  (§1)  can  naturally  produce  a  cooler 
infrared  color,  because  the  contribution  from  an  outer,  cooler  com¬ 
ponent  becomes  important  and  makes  the  optical  Seyfert  signa¬ 
tures  fainter  or  even  vanish,  because  even  the  torus  axis  direction 
can  be  opaque  to  the  bulk  of  the  AGN’s  ionizing  radiation. 

5.5.  Future  Prospects 

Although  infrared  spectroscopic  energy  diagnostics  enable 
us  to  discuss  the  overall  properties  of  ULIRGs  statistically, 
follow-up  observations  of  infrared-selected  buried  AGN  candi¬ 
dates  are  very  important  to  cross-check  the  reliability  of  these 
methods.  For  example,  we  argue  an  exceptionally  centrally  con¬ 
centrated  extreme  starburst  is  unlikely,  based  on  the  emission 
surface  brightness  and  absolute  luminosity  (§  5.2.5).  X-ray  obser¬ 
vations  can  be  particularly  powerful  in  distinguishing  a  buried 
AGN  from  an  extreme  starburst,  because  a  luminous  AGN  is  a  much 
stronger  X-ray  emitter  than  any  kind  of  starburst.  Since  the  column 
densities  in  the  cores  of  ULIRGs  can  be  large  enough  to  be  Comp¬ 
ton  thick  (X-ray  absorption  Ah  >  1024  cm-2),  high-sensitivity 
X-ray  observations  at  E  >  10  keV  are  particularly  desirable. 


98 


IMANISH1  ET  AL. 


Vol.  171 


In  a  strongly  X-ray-emitting  buried  AGN  surrounded  by  dense 
gas  and  dust,  X-ray  dissociation  regions  (XDRs;  Maloney  et  al. 
1996)  should  develop,  rather  than  the  photodissociation  regions 
( PDRs)  usually  seen  in  a  starburst  (a  strong  UV emitter).  Finding 
XDR  signatures  can  thus  be  another  powerful  way  to  test  the 
buried  AGN  scenario.  Imanishi  et  al.  (2004,  2006b)  found  that 
several  infrared-selected  buried  AGN  candidates  tend  to  show 
millimeter  molecular  line  flux  ratios  expected  from  XDRs.  For 
one  of  the  strongest  infrared-selected  buried  AGN  candidates, 
IRAS  08572+3915,  Shirahata  et  al.  (2007)  detected  CO  (v  =  1-0) 
absorption  features  up  to  J  =  17  in  ground-based  infrared  4.5- 
5  /im  observations.  It  is  argued  that  such  high-7  lines  can  be  ex¬ 
plained  only  by  XDRs  (Spaans  2006),  further  supporting  the 
buried  AGN  scenario.  A  further  search  for  XDR  signatures  in 
ULIRGs  may  reinforce  the  buried  AGN  picture  suggested  from 
the  infrared  spectroscopy. 

6.  SUMMARY 

We  presented  Spitzer  IRS  5-35  pm  low-resolution  (R  ~  100) 
spectra  of  ULIRGs  at  z  <  0.15  in  the  IRAS  1  Jy  sample.  ULIRGs 
classified  optically  as  non-Seyferts  (LINERs  and  FI  n  regions), 
which  comprise  70%  of  the  1  Jy  sample  at  z  <  0.15,  were  our 
prime  targets  to  search  for  signatures  of  optically  undetected 
AGNs  deeply  buried  in  dust  along  virtually  all  lines  of  sight, 
without  well-developed  narrow-line  regions.  In  total,  28  LINER 
and  20  H  n  region  ULIRGs  (a  complete  sample)  were  observed. 
Three  methods  were  mainly  utilized:  (1)  the  equivalent  widths  of 
the  PAH  emission  features  at  6.2,  7.7,  and  11.3  /tm,  (2)  the  ab¬ 
solute  optical  depth  of  the  9.7  /tm  silicate  dust  absorption  feature, 
and  (3)  the  9.7-18  /im  silicate  dust  absorption  optical  depth  ra¬ 
tio.  The  first  method  traces  the  hardness  of  the  radiation  field 
(starburst-like  soft  or  AGN-like  hard).  The  latter  two  methods 
probe  the  geometry,  i.e.,  whether  the  dust  and  energy  sources 
are  spatially  well  mixed  (a  normal  starburst)  or  the  energy  source 
is  more  centrally  concentrated  than  dust  (a  buried  AGN).  We 
arrived  at  the  following  main  conclusions. 

1 .  For  the  majority  of  ULIRGs,  PAH  emission  features  were 
detected,  suggesting  the  presence  of  PAH-emitting  normal  star- 
burst  activity.  However,  the  observed  PAH-to-infrared  luminosity 
ratios  were  systematically  smaller  than  in  modestly  obscured 
(Ay  <  20-30  mag)  starburst  galaxies.  Taken  at  face  value,  the 
PAH-emitting  starbursts  with  modest  dust  obscuration  are  ener¬ 
getically  insignificant  in  ULIRGs.  The  dominant  energy  sources 
are  (1)  very  highly  obscured  (Ay  ys>  30  mag)  PAH-emitting  star- 
bursts  with  mixed  dust /source  geometry,  whose  PAH  flux  is  at¬ 
tenuated,  and/or  (2)  non-PAH-emitting  AGNs. 

2.  In  19  out  of  the  48  (40%)  ULIRGs  we  studied,  the  ob¬ 
served  PAH  equivalent  widths  were  found  to  be  substantially 
lower  than  normal  starburst  galaxies  at  either  of  the  6.2,  7.7,  or 
11.3  pm  PAH  features.  This  fraction  increased  to  22/48  (46%)  if 
we  adopted  the  6.2  p m  PAH  equivalent  widths  after  correction 
for  the  ice  and  HAC  absorption  features.  The  highly  obscured 
energy  sources  of  these  ULIRGs’  cores  were  argued  to  be  domi¬ 
nated  by  non-PAH-emitting  AGNs,  rather  than  highly  obscured 
PAH-emitting  starbursts. 

3.  In  20/48  (42%)  of  ULIRGs,  the  observed  optical  depths 
of  the  9.7  /im  silicate  dust  absorption  feature  (Tg7)  were  sub¬ 
stantially  larger  than  the  upper  limit  determined  by  mixed  dust/ 
source  geometry  in  normal  starbursts.  A  non-PAH-emitting  buried 
AGN  with  centrally  concentrated  energy  source  geometry  is  the 


most  natural  explanation  for  the  energy  sources  at  the  cores  of 
these  ULIRGs.  Comparisons  of  the  r)  7  value  with  the  PAH 
equivalent  widths,  and  with  the  11.3-6.2  pm  PAH  luminosity 
ratio,  supported  the  buried  AGN  scenario  in  these  ULIRGs. 

4.  For  selected  ULIRGs  with  small  PAH  contamination 
(mostly  ULIRGs  with  large  Tg  7),  the  18-9.7  pm  silicate  absorp¬ 
tion  optical  depth  ratios  were  often  significantly  smaller  than  the 
intrinsic  ratio.  A  strong  dust  temperature  gradient  around  a  cen¬ 
trally  concentrated  energy  source  geometry  was  strongly  sug¬ 
gested,  further  supporting  the  buried  AGN  scenario. 

5.  Sixteen  out  of  48  (33%)  of  optical  non-Seyfert  ULIRGs 
showed  strong  signatures  of  luminous  buried  AGNs,  and  26/48 
(54%)  ULIRGs  displayed  some  signatures.  Given  that  30%  of 
ULIRGs  showed  optical  Seyfert  signatures  and  were  known  to 
possess  luminous  AGNs  surrounded  by  dusty  tori,  it  was  sug¬ 
gested  that  luminous  AGNs  are  present  in  at  least  50%,  and  pos¬ 
sibly  >65%,  of  the  complete  ULIRG  sample  at  z  <  0.15. 

6.  The  fraction  of  ULIRGs  showing  luminous  buried  AGN 
signatures  was  always  higher  in  optical  LINER  ULIRGs  than  H  n 
region  ULIRGs  using  any  of  the  above  methods.  This  was  quali¬ 
tatively  explained  by  the  scenario  that  the  starburst  regions  probed 
in  the  optical  are  dustier  in  LINER  ULIRGs  than  in  H  n  region 
ULIRGs. 

7.  For  ULIRGs  with  strong  buried  AGN  signatures  and  small 
PAH  contamination,  the  absorption-corrected  intrinsic  luminos¬ 
ities  of  the  buried  AGNs  were  found  to  be  quantitatively  com¬ 
parable  to  1012  Lq,  accounting  for  the  bulk  of  the  observed 
infrared  luminosities  of  ULIRGs. 

8.  Most  of  the  luminous  buried  AGN  candidates  were  found 
in  ULIRGs  with  cool  far-infrared  colors. 

9.  The  above  conclusions  drawn  from  the  Spitzer  IRS  5— 
35  /tm  spectra  were  generally  consistent  with  those  based  on 
3-4  /tm  spectra  obtained  by  Imanishi  et  al.  (2006a).  Our  main 
conclusions  were  that  optically  undetected  luminous  buried  AGNs 
are  surely  present  in  a  significant  fraction  of  non-Seyfert  ULIRGs. 
It  is  therefore  very  important  to  carefully  and  quantitatively  in¬ 
vestigate  the  energetic  roles  of  optically  elusive  luminous  bur¬ 
ied  AGN  s,  if  we  are  to  understand  the  true  nature  of  the  nearby 
ULIRG  population. 
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APPENDIX 

THE  MAXIMUM  t/7  VALUE  IN  A  MIXED  DUST/SOURCE  GEOMETRY 
Figure  1  a  shows  a  schematic  diagram  of  a  typical  mixed  dust/source  geometry.  In  this  geometry,  the  observed  flux  7(r„)  is  given  by 

i(n)  -  70(a) x1-^,  (Al) 

T"A 


where  Iq(X)  is  an  unattenuated  intrinsic  flux  and  is  the  optical  depth  at  each  wavelength,  which  takes  different  values  inside  ab¬ 
sorption  features  and  nearby  continua  outside  the  features. 

We  here  consider  radiation  with  an  intrinsically  smooth  spectral  energy  distribution,  namely  7o(8  /t m)  ~  /o(9.7  /tm)  ~  /o(13  /tm), 
such  as  expected  for  dust  emission  in  ULIRGs.  The  peak  absorption  optical  depth  of  the  9.7  /tm  silicate  feature,  relative  to  the  nearby 
continuum  at  8  and  13  /tm  (which  corresponds  to  r9  7  and  not  r9  7  in  our  definition  in  §  4.2),  is 


^  J  _  q  ^cont 

7"97sil 

V  Tcont 

\  —  T"97sil 
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1  —  e_T97sil 

Tcont  / 
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where  t97s;i  and  rcont  refer  to  the  optical  depths  at  the  peak  of  the  9.7  /im  silicate  absorption  feature  and  nearby  continuum  at  8  and 
13  /tm  just  outside  the  feature,  respectively.  The  t97s;i  value  is  always  larger  than  the  rcont  value,  because  r97sii  is  the  superposition 
of  the  optical  depth  for  continuum  emission  and  the  9.7  /im  silicate  feature. 

In  the  case  of  the  Galactic  dust  extinction  curve  derived  by  Rieke  &  Lebofsky  (1985), 

Tcont  =  0. 025/1. 08Hf/,  (A4) 

T9is\\  =  0.087 / 1 .08Ay,  (A5) 


where,  by  definition,  Ay  =  1.08  x  Ty. 

The  absorption  strength  of  the  9.7  /tm  silicate  feature,  relative  to  Ay,  is  different  by  a  factor  of  ~2  in  different  directions  in  the 
Galaxy  (Roche  &  Aitken  1984,  1985).  If  we  take  this  into  account, 

Tcont  =  0.025/ 1 ,08Ay,  (A6) 

r97sil  =  (0.05-0.  l)Ay  +  0.025/1.08 Ay.  (A7) 

The  rcont  and  t97s;i  values  can  be  arbitrarily  large  with  increasing  amount  of  dust  (Ay)  in  the  mixed  dust/source  geometry.  The  r9  7 
value  is  an  increasing  function  of  Ay,  and  the  first  term  in  the  equation  (A3)  [(1  —  e-r“n,)/(l  —  e~T97si1)]  becomes  unity  for  a  very 
large  Ay  value.  Thus,  the  upper  limit  to  r97  is 

Tg  7(upper  limit)  =  In  j  T97sl1  )  _  (A8) 

\  ^"cont  / 


We  obtain  a  stringent  upper  limit  of  r97  ~  1.7  in  this  mixed  dust/source  geometry  for  any  amount  of  dust. 

The  TCont/A  v  and  r97siiM  y  ratios  can  vary  substantially  with  varying  dust  grain  size,  because  the  wavelength  difference  between  the 
optical  Kband  and  infrared  ~  1 0  /tm  is  large.  However,  the  upper  limit  of  r9  7  in  the  mixed  dust/source  geometry  is  dependent  only  on 
Tcont /T97sii,  which  should  be  less  sensitive  to  varying  dust  grain  size  (see  Fig.  6b),  because  their  wavelengths  are  very  close.  Therefore, 
the  above  conclusion  about  the  upper  limit  of  r9  7  in  the  mixed  dust/source  geometry  should  be  robust  to  possible  differences  in  dust 
grain  size  between  the  cores  of  ULIRGs  and  the  Galactic  diffuse  interstellar  medium. 
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